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Let ({^i(t)}igz<i)t>o be the system of interacting diffusions on 
[0, oo) defined by the following collection of coupled stochastic differ- 
ential equations: 

dX,{t) = a{i,j)[Xj{t) - X,{t)]dt+ y^bX,{t)^ dW,{t), 

i£Z'^,t>0. 

Here, •) is an irreducible random walk transition kernel on Z'^ x 
Z'', b£ {0,oo) is a diffusion parameter, and {{Wi{t)}i^^d)t>o is a 
collection of independent standard Brownian motions on R. The ini- 
tial condition is chosen such that {Xi(0)}i^^d is a shift-invariant and 
shift-ergodic random field on [0, oo) with mean O £ (0, oo) (the evo- 
lution preserves the mean). 

We show that the long-time behavior of this system is the re- 
sult of a delicate interplay between a(-,-) and b, in contrast to sys- 
tems where the diffusion function is subquadratic. In particular, let 
^ihj) ~ ^[o-ihj) + o-Uji-)]; ijj £ ^'') denote the symmetrized transi- 
tion kernel. We show that: 

(A) If a(-, •) is recurrent, then for any b > the system locally 
dies out. 

(B) If a(-, •) is transient, then there exist fe* > 62 > such that: 
(Bl) The system converges to an equilibrium (with mean 

e) ff 0<fo<fo,. 

(B2) The system locally dies out if 6 > 6, . 

(B3) has a finite 2nd moment if and only if < b < 62. 
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(B4) The 2nd moment diverges exponentially fast if and 
only if 6 > &2 • 

The equilibrium i^e is shown to be associated and mixing for all 
< 6 < 6*. We argue in favor of the conjecture that 6, > fe2- We 
further conjecture that the system locally dies out at 6 = 6,. 

For the case where a(-,-) is symmetric and transient we further 
show that: 

(C) There exists a sequence 62 > 63 > 64 > ■ ■ ■ > such that: 

(CI) U0 has a finite mth moment if and only if < 6 < 6m- 
(C2) The mth moment diverges exponentially fast if and 
only if 6 > 6m. 

(C3) 62 < (m - l)bm < 2. 

(C4) limm^oo (m — l)6m = c = sup^>2 irn—l)bm- 
The proof of these results is based on self-duality and on a represen- 
tation formula through which the moments of the components are 
related to exponential moments of the collision local time of random 
walks. Via large deviation theory, the latter lead to variational ex- 
pressions for 6, and the 6m 's, from which sharp bounds are deduced. 
The critical value 6, arises from a stochastic representation of the 
Palm distribution of the system. 

The special case where a(-, •) is simple random walk is commonly 
referred to as the parabolic Anderson model with Brownian noise. 
This case was studied in the memoir by Carmona and Molchanov 
[Parabolic Anderson Problem and Intermittency (1994) Amer. Math. 
Soc, Providence, RI], where part of our results were already estab- 
lished. 

1. Introduction and main results. 

1.1. Motivation and background. This paper is concerned with the long- 
time behavior of a particular class of systems with interacting components. 
In this class, the components are interacting diffusions that take values in 
[0, 00) and that are labelled by a countably infinite Abelian group /. The 
reason for studying these systems is two- fold: new phenomena occur, and a 
number of methodological problems can be tackled that are unresolved in 
the broader context of interacting systems with noncompact components. 
We begin by describing in more detail the background of the questions to 
be addressed. 

A large class of interacting systems has the property that single compo- 
nents change according to a certain random evolution, while the interaction 
between the components is linear and can be interpreted as migration of 
mass, charge or particles. Examples are: 

(1) Interacting particle systems, for example, voter model [34], branching 
random walk [22, 36], generalized potlatch and smoothing process [35], 
binary path process [31], coupled branching process [28, 29], locally de- 
pendent branching process [3], catalytic branching [27, 37]. 
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(2) Interacting dijfusions, for example, Fisher-Wright diffusion [9, 10, 13, 
14, 24, 25, 32, 33, 40, 41, 44], critical Ornstein-Uhlenbeck process [19, 
20], Feller's branching diffusion [15, 41], parabolic Anderson model with 
Brownian noise [7]. 

(3) Interacting measure-valued diffusions, for example, Fleming-Viot pro- 
cess [17], mutually catalytic diffusions [18], catalytic interacting diffu- 
sions [30]. 

Most of these systems display the following universality: independently of 
the nature of the random evolution of single components, the ergodic be- 
havior of the system depends only on recurrence versus transience of the 
migration mechanism. More precisely, if the symmetrized migration kernel 
is recurrent then the system approaches trivial equilibria (concentrated on 
the "traps" of the system) , whereas if the symmetrized migration kernel is 
transient then nontrivial extremal equilibria exist that can be parametrized 
by the spatial density of the components. 

In this paper we study an example in a different universality class, one 
where the nature of the random evolution of single components does influ- 
ence in a crucial way the long-time behavior of the system. In particular, 
we consider a system where the components evolve as diffusions on [0, oo) 
with diffusion function bx'^ and interact linearly according to a random walk 
transition kernel. Such a system is called the parabolic Anderson model with 
Brownian noise in the special case where the random walk is simple. In the 
recurrent case the system, as before, approaches a trivial equilibrium (con- 
centrated on the "trap" with all components 0), so local extinction prevails. 
However, in the transient case we find three regimes, separated by critical 
thresholds 6* > 62 > (see Figure 1): 

(I) ("low noise") < 5 < 62: equilibria with finite 2nd moment. 

(II) ( "moderate noise" ) 62 < 6 < 6* : equilibria with finite 1st moment and 
infinite 2nd moment. 

(III) ("high noise") 6 > 6*: local extinction. 

We will show that the strict inequality 6* > 62 depends on a large deviation 
principle for a renewal process in a random environment. This large deviation 
principle will be addressed in a forthcoming paper (Birkner, Greven and 
Hollander [5]). Local extinction at & = is a subtle issue that remains open. 

For the case where the random walk transition kernel is symmetric we do 
a finer analysis. We show that in regime (I) there exists a sequence 62 ^ ^3 ^ 
64 > • • • > such that the equilibria have a finite mth moment if and only if 
<b <bm, while the mth moment diverges exponentially fast if and only if 
b> bm (see Figure 1). Moreover, we show that 62 < (nT- — l)^m < 2 and that 
limm->oo("T' — l)^m = c = sup^>2(m — l)6m- We show that in regimes (I) and 
(II) the equilibria are associated and mixing. We show that the critical value 
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6=1, separating regimes (II) and (III) is linked to the Palm distribution of the 
system. 

The reason for the above phase diagram is that there are two compet- 
ing mechanisms: the migration pushes the components toward the mean 
value of the initial configuration, while the diffusion pushes them toward 
the boundary of the state space. Hence, there is a dichotomy in that either 
the migration dominates (giving nontrivial equilibria) or the diffusion dom- 
inates (giving local extinction). In the class of interacting diffusions we are 
concerned with here, the migration and the diffusion have a strength of the 
same order of magnitude and therefore the precise value of the diffusion pa- 
rameter in relation to the migration kernel is crucial for the ergodic behavior 
of the system. 

Our results are a completion and a generalization of the results in the 
memoir of Carmona and Molchanov [7]. In [7], Chapter III, the focus is on 
the annealed Lyapunov exponents for simple random walk, that is, on Xm{b), 
the exponential growth rate of the mth moment of Xo(t), for successive m. 
It is shown that for each m there is a critical value bm where Xm{b) changes 
from being zero to being positive (see Figure 2), and that the sequence {bm) 
has the qualitative properties mentioned earlier, that is, 6^ = for all m in 
d= 1,2 (recurrent case) and 62 ^ 63 > 64 > • • • > in d > 3 (transient case). 
No existence of and convergence to equilibria is established below b2, nor 
is any information on the equilibria obtained. There is also no analysis of 
what happens at the critical values. In our paper we are able to handle these 
issues due to the fact that we have variational expressions for Xm{b) and bm, 
which give us better control. In addition, we are able to get sharp bounds 
on bm that are valid for arbitrary symmetric random walk, which results in 
strict inequalities between the first few 6m 's- 

In [7], Chapter IV, an analysis is given of the quenched Lyapunov exponent 
for simple random walk, that is, on x*ib), the a.s. exponential growth rate 
of Xo^t). It is shown that x*ib) is negative for all 6 > in d = 1, 2 (recurrent 
case) , negative for b>b^, and zero for < 6 < 6* in d > 3 (transient case) for 
some 6^, > 62 (see Figure 2). This corresponds to the crossover at 6^,, except 
for the proof that 6^, > 62 , which we defer to a forthcoming paper [5] . In [7] , 
Chapter IV, it is further shown that x* (b) has a singular asymptotics for b — > 
00. This asymptotics has been sharpened in a sequence of subsequent papers 



(I) 



(11) 



(III) 



Fig. 1. Phase diagram for the transient case. 
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by Carmona, Molchanov and Viens [8], Carmona, Koralov and Molchanov 
[6] and Cranston, Mountford and Shiga [12]. 

A scenario as described above is expected to hold for a number of interact- 
ing systems where the components take values in a noncompact state space, 
for example, generalized potlatch and smoothing [35] and coupled branching 
[28, 29]. But for none of these systems has the scenario actually been fully 
proven. 

1.2. Open problems. We formulate a number of open problems that are 
not addressed in the present paper: 

(A) Show that 62 > ^3 > ^4 > • • • ■ This property is claimed in [7], Chapter 
III, Section 1.6, but no proof is provided. We are able to show that 
b2> bs> ■ ■ ■ > bm for an arbitrary symmetric random walk for which 
the average number of returns to the origin is < l/(m — 2). For m = 3, 
this includes simple random walk in d>3. 

(B) Show that the system locally dies out at the critical value b^. 

(C) Show that x*ib) < for 6 > 6^,, that is, show that there is no intermedi- 
ate regime where the system locally dies out but only subexponentially 
fast. Shiga [41] has shown that the system locally dies out exponentially 
fast for b sufficiently large. 

(D) Find out whether there exists a characterization of 6^, in terms of the 
collision local time of random walks. This turns out to be a subtle 
problem, which has analogues in other models (see Birkner [3]). We 
find that such a characterization does exist for bm and for a certain 
with bif > 6=n<. We have a characterization of 6* in terms of the Palm 
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distribution of our process, but this is relatively inaccessible. It therefore 
is a subtle problem to decide whether 6^, = 6** or 6* > 6** . 

1.3. Outline. In Section 1.4 we define the model, formulate a theorem by 
Shiga and Shimizu [42] stating that our system of interacting diffusions has 
a unique strong solution with the Feller property, and introduce some key 
notions. In Section 1.5 we formulate two more theorems, due to Shiga [41] 
and to Cox and Greven [10], stating that our system locally dies out in the 
recurrent case and has associated mixing equilibria with finite 2nd moment 
in the transient case in regime (I). We complement these two theorems with 
two new results, stating that our system has associated mixing equilibria 
with finite 1st moment in the transient case in regime (II) and no equilibria 
in the transient case in regime (III). In Section 1.6 we present our finer 
results for regime (I), and have a closer look at regimes (II) and (III) as 
well, although much less detailed information is obtained for these regimes. 

Sections 2-4 contain the proofs. Section 2 is devoted to moment calcula- 
tions, which are based on a (Feynman-Kac type) representation formula for 
the solution of our system due to Shiga [41]. Through this representation 
formula, we express the moments of the components of our system in terms 
of exponential moments of the collision local time of random walks. Through 
the latter we are able to establish convergence to a (possibly trivial) equi- 
librium and to prove that this equilibrium is shift invariant, ergodic and 
associated. In Section 3 we study the exponential moments of the collision 
local time with the help of large deviation theory, which leads to a detailed 
analysis of the critical thresholds b function of m in regime (I) , as well 

as to a description of the behavior of the system at hm- Section 4 looks at sur- 
vival versus extinction and relates the critical threshold h-f between regimes 
(II) and (III) to the so-called Palm distribution of our system, where the 
law of the process is changed by size biasing with the value of the coordinate 
at the origin. There we argue in favor of the strict inequality 6^, > 62 1 which 
relies on an explicit representation formula for the Palm distribution. 

1.4. The model. The models that we consider are systems of interacting 
diffusions X = (X{t))t>o, where 

(1.1) X{t)={X,{t)}ieie[0,^Y, 

with / a countably infinite Abelian group. The evolution is defined by the 
following system of stochastic differential equations (SSDE): 

dXi{t) = J2 a{i,j)[Xj{t) - X,{t)] dt + ^JbXi{t)^ dW^{t), 

(1.2) 

i£l,t>0. 

Here: 
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(i) a{-, •) is a Markov transition kernel on / x /. 

(ii) b G (0, oo) is a parameter. 

(iii) W = {{Wi{t)}i^[)t>o is a collection of independent standard Brown- 
ian motions on M. 

Equation (1.2) arises as the continuum limit of a self-catalyzing branching 
Markov chain whose branching rate depends on the local population size. 
As initial condition we take 

(1.3) x(o)g£:i, 

where 

(1.4) = jrr = (xjie/ G [0, 00)^^7^^^ < ooj C ^^(7) 
for any 7 = (7j)jg/ satisfying the requirements 

7/ > yiei, 

(1.5) ^7i<oo, 
is/ 

3M < cx) : ^ 7ia{i,j) < M-fj Vj G /. 

We endow <Si with the product topology of [0, 00)^ (see Liggett and Spitzer 
[38]). 

Since |/| = cx), it is not possible to define the process uniquely in the strong 
sense on [0, 00)^ without putting growth conditions on the initial configura- 
tion, as in (1.4). However, the dependence of £1 on 7 is not very serious. For 
example, every probability measure p on [0, 00)^ satisfying sup^^j EP{Xi) < 
00 is concentrated on £1 regardless of the 7 chosen {E^ denotes expectation 
with respect to p). We also need the space £2 C L'^i'y), which is defined as in 
(1.4) but with the condition J^ieili^i < cxd replaced by lii^i)'^ < 00. 

The most basic facts about the process {X{t))t>o are summarized in the 
following result. 

Theorem 1.1 (Shiga and Shimizu [42]). 

(a) The SSDE in (1.2) has a unique strong solution {X{t))t>o on £1 
with continuous paths. 

(b) {X{t))t>o is the unique Markov process on£i whose semigroup {S{t))t>o 
satisfies 

(1.6) S{t)f-f= ['s{s)Lfds, feCl{£i), 

Jo 
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where Cq{£i) is the space of functions on £i depending on finitely many 
components and twice continuously differentiahle in each component, and L 
is the pregenerator 

{\ df 1 9^ / 

(1.7) 

(c) Restricted to £2, {X{t))t>o is a diffusion process with the Feller prop- 
erty. 

The model defined by (1.2) represents a special case of the SSDE 
dX,{t) = Y,a{i,j)[Xj{t) - X,{t)] dt + ^jg{X,{t)) dWi{t), 

(1.8) 

iei,t>o, 

with g: (-co, 00) — > [0,cx)) some locally Lipschitz continuous function. This 
SSDE has, as far as its long-time behavior is concerned, four important 
classes: 

(i) g{x)>Oon (0,1). 

Examples: g{x) = x{l — x) Fisher- Wright, g{x) = {x(l — x))^ Ohta-Kimura. 

(ii) g{x) > on (—00, 00) and g{x) = o(x^) as x — > ±00. 
Example: ^(x) = a'^ critical Ornstein-Uhlenbeck. 

(iii) g{x) > on (0, 00) and ^(x) = o(x^) as x — > 00. 
Example: ^(x) = x Feller's continuous-state branching diffusion. 

(iv) g{x) > on (0, 00) and g{x) ~ 6x^ as x — > 00. 

Classes (i)~(iii) are well understood [9, 10, 19, 20, 24, 40, 41]. The qualitative 
properties of the process defined by (1.8) are similar for these three classes, 
and the universality of the long-time behavior as a function of g has been 
systematically investigated via renormalization methods [1, 2, 13, 14, 15, 16, 
26]. Class (iv), which is the subject of the current paper, is very different. 
For the case where a(-, •) is simple random walk, this class was investigated 
in [41] and in the memoir by Carmona and Molchanov [7], where some of 
our results were already established. 

The long-time behavior of the process defined by (1.2) is fairly complex. 
In order to keep the exposition transparent, we restrict our analysis to a 
subclass of models given by the following additional requirements: 

I = Z'^, d>l, 

(1.9) a(-, •) is homogeneous: a{i,j) = a(0, j — i) G /, 

00 

a(-, •) is irreducible: ''^[a"{i,j)-\-a"{j,i)]>0 \/i,jGl. 

n=0 
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Moreover, we put a(0, 0) = 0. 

Before we start, let us fix some notation. We write V{£i) for the set of 
probability measures on {£i,13{£i)), with B the Borel u-algebra. For p£ 
V{£i), we write to denote expectation with respect to p. A measure 
yO G 'Pi£i) is called shift-invariant if 

(1.10) p((X,).e7G^)=/>((X,+,■)^6/GA) yj£iyAeB{£i), 
is called mixing if 

(1.11) hm EP{f[goak]) = EP{f)E''{g) 

for all bounded f,g:£i that depend on finitely many coordinates, where 
cTfc is the A;-shift acting on /, and is called associated if 

(1.12) EP{hh)>EP{h)EP{h) 

for all bounded /i, /2 : <?i ^ M that depend on finitely many coordinates and 
that are nondecreasing in each coordinate. 
We further need 

T = {p G ^{£1) : p is shift-invariant}, 

(1.13) 

T^ = {p(^T:EP{X^)< 00}, 

and 

(1.14) = {pGT^: pis shift ergodic,^^(Xo) = e}, eG[0,oo). 

The set of extreme points of a convex set C is written Cg. The element 
{xi)i^i with Xi = Q for alH € / is denoted by G. The initial distribution of 
our system is denoted hy p = £(X(0)) and is assumed to be concentrated 
on £1 . The symbols P, E without index denote probability and expectation 
with respect to p and the Brownian motion driving (1.2). The notation w-lim 
means weak limit. 



1.5. Phase transitions. In Theorems 1.2-1.5 below we state our main 
results on the long-time behavior of {X{t))t>o and on the properties of its 
equilibria. Let 

(1.15) X = {p G 'P((5i) : p is invariant} 

be the set of all equilibrium measures p of (1.2), that is, pS{t) = p for all 
t > 0. This set of course depends on a(-, •) and b. 
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1.5.1. Recurrent case. The ergodic behavior of our system is simple when 
a(-, •) defined by 

(1.16) a{i,j) = ^[a{i,j)+a{j,i)], 

is recurrent. Namely, the process becomes extinct independently of the value 
of b. 

Theorem 1.2 (Shiga [41]). Ifa{-,-) is recurrent, then for every b>0 
and every initial distribution : 

(1.17) w-\imC{X{t)) = 6o. 

Consequently, there exists no equilibrium in other than Sq, that is, 

(1.18) J nTi = Jo- 
Using the fact that if /i G Tq then E{Xi{t)) = @ for alH G / and t > 0, 

we conclude from Theorem 1.2 that the system clusters, that is, on only few 
sites there is a nontrivial mass but at these sites the mass is very large (for 
t large). 

1.5.2. Transient case: regimes (I), (II) and (III). In the case where a (•, •) 
is transient, the ergodic behavior of our system depends on the parameter 
b and we observe interesting phase transitions. There are three regimes, 
separated by two critical values. 

(I) Small b. Define the Green function 

oo 

(1.19) G{i,j) = J2a^{i,j), i,j€l, 

n=0 

and put 

(1.20) h 



G(0,0) 

We first consider the regime 

(1.21) (I) a(-,-) transient, b£{QM)- 

Theorem 1.3 (Shiga [41], Cox and Greven [10]). In regime (I): 
(a) For fi = 6q with Q G [0,oo) the following limit exists: 

(1.22) 1/0 =w-lim£(X(t)). 
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(b) The measure satisfies 

vq is shift- invariant, mixing and associated, 

(1.23) 

(Xq) = 0; z^e is not a point mass if @ > 0, 
E''^{[Xo]^)<oo. 

(c) The set of shift- invariant extremal equilibria is given by 
(1-24) (XnTl)e = {z^0}e6[o,oo). 

(d) For every fi £Tq with Q G [0, oo): 

(1.25) w-lim£(X(i)) = i/e. 

t — >oo 

(e) For every fi with E^{Xq) = cxo; 

(1.26) w-lim£(X(t)) = 5oo. 

t— >oo — 

Consequently, 

(1.27) {lr^T)e = {lr^T^)e. 

Theorem 1.3 tells us that if b remains below an a(-, •)-dependent threshold, 
then the process {X[t))t>Q exhibits persistent behavior, in the sense that an 
equilibrium is approached with a spatial density equal to the initial spatial 
density and with a one-dimensional marginal that has a finite 2nd moment. 
This equilibrium is nontrivial unless the initial state is identically 0. If, on 
the other hand, the initial spatial density is infinite, then every component 
diverges in probability. 

(II) Moderate b. We next consider the regime 

(1.28) (II) a(-,-) transient, 6 £[62,6*). 

In Section 4 we will obtain a variational expression for 6* [see (4.19)]. This 
expression will turn out to be somewhat delicate to analyze. 

Theorem 1.4. In regime (II); 

(a) The same properties hold as in Theorem 1.3(a) and (c)-(e). 

(b) The measure vq satisfies 

i/Q is shift-invariant, mixing and associated, 

(1.29) 

(Xq) =@;uq is not a point mass if @ > 0, 
E'^^{[Xo]^) = ^. 
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Theorem 1.4, which will be proved in Section 2, says that for moderate 
b the equilibria continue to exist and to be well behaved, but with a one- 
dimensional marginal having infinite 2nd moment. The latter has important 
consequences for the fluctuations of the equilibrium in large blocks. Indeed, 
in regime (I) we may expect Gaussian limits after suitable scaling (see, e.g., 
Zahle [46, 47] in a different context), while in regime (II) we may expect 
non-Gaussian limits. In regime (II), the tail of Xq under is likely to be 
of stable law type, but a closer investigation of this question is beyond the 
scope of the present paper. 

(Ill) Large b. Finally, we consider the regime 

(1.30) (III) a (•,•) transient, b£[b^,oo). 

Theorem 1.5. In the interior of regime (III), for every : 

(1.31) w-limC{X{t)) = 6o. 

t—fOC ~ 

Consequently, I HT^ = Sq. 

We conjecture that there is local extinction at b = b^. 

Theorem 1.5, which will be proven in Section 4, shows that for large b 
again clustering occurs, that is, the same situation as described in Theorem 
1.2 for the case where a(-,-) is recurrent. 

1.6. Finer analysis of the transient case. In Section 1.5 we saw that 
different values of b lead to qualitatively different behavior of the process 
{X{t))t>Q- Therefore the question arises in which way the value of b influ- 
ences the properties of the process within one regime. For part of this finer 
analysis we need to assume that a(-,-) is symmetric: 

(1.32) a{i,j) = a{j,i) yijel. 

1.6.1. Regime (I). Let = {(,{t))t>o be the random walk on / with tran- 
sition kernel a(-,-) and jump rate 1, starting at 0. For m > 2, let = 
(?i) • • • )Cm) be m independent copies of ^, and define the differences random 
walk ry(™) = (t?*-™^ (t))t>o by putting 

(1-33) ^^™^(i) = (ep(t)-e,(t))i<p<,<™. 

This is a random walk on the subgroup of /(V2)m(m-i) generated by 

all the possible pairwise differences of m elements of I, with jump rate m 
and transition kernel a^"^\-, •) that can be formally written out as 

aM(x,y) = a(")(0,y-x) 

(1.34) 



-1 lib 

-Y^l{jDr=y-x} 

m , 



x,y 
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where Dr is the triangular array of —1,0, +l's given by 

(1.35) Dr = {6pr — Sqr)l<p<^q<m 

and jDr denotes the triangular array obtained from D,. by multiplying all 
its elements with the vector j . The factor ^ comes from the fact that the m 
random walks jump one at a time. Note that a^"^^(-, •) is symmetric because 
of our assumption in (1.32). Note that a^'^\-,-) =a(-,-), the symmetrized 
transition kernel defined in (1.16), which is symmetric even without (1.32). 
The differences random walk is to be seen as the evolution of the random 
walks "relative to their center of mass." This will serve us later on. 
Define the Green function 

CXj 

(1.36) G^"'\x,y) = ^[a(-)]"(x,y), x,yG 

n=0 

Also define the collision function (J^™) : J^™) Nq as 
(1.37) 



l<p<q<m 

■2- — {^p ■2-g)l<p<(jr<m; ^p, £ 

and put 

(1.38) S('") = supp(tt(™))c 

Define 



(1.39) K("'\x,y) = ,fi^)G^'^H^,y)\/¥^), x,y&S^^\ 

Viewed as an operator acting on £'^{S^'^^), K^"^^ is self-adjoint, positive and 
bounded. The latter two properties will be proved in Section 2. 

The following result shows that in regime (I) there is an infinite sequence 
of critical values characterizing the convergence of successive moments. 

Theorem 1.6. Suppose that a(-,-) is symmetric. Then, in regime (I), 
there exists a sequence 62 > 63 > &4 > • • • such that: 

(a) If = 6@ with Q G (0, 00), then 

< 00, for b < bm, 



t-»oo = 00, forb>b. 

(b) // /_f = (50 with G (0, cxd), then 

(1.41) hm - \ogE{[Xo{t)r) = Xm{b) 

t~>oo t 

exists with 

0, forb<bm, 



(1-42) Xm{b) 



> 0, for b>b, 



m • 
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(c) The critical value bm has the representation 

in 

(1.43) bm = ^ 

Xm 

with Am G (0, oo) the spectral radius of K^"^^ in i'^{S^"^^). This spectral radius 
is an eigenvalue if and only if bm-i > bm- 

(d) The critical value 62 is given by (1.20), and 

(1.44) 62 > 63 > ^4 > • • • > 0. 

Moreover, 

(1-45) =b,<{m- l)bm < < 2, 

and limm^oo(?^ — l)^m exists. 

(e) The function b^ m^"^^^^ ^"^ convex on [0,oo) and strictly increasing 
on [bm,oo), with 

(1.46) Urn ^xmib) = lim-l). 



Theorem 1.6 will be proven in Section 3. Part (a) tells us that equilibria 
with finite mth moment exist if and only ii <b <bm- Part (b) tells us that 
the mth moment diverges exponentially fast if and only if 6 > 6^- The limit 
Xmib) is the mth annealed Lyapunov exponent. Part (c) gives a variational 
representation for bm- Part (d) gives sharp bounds for bm and shows that 
the tail of the one-dimensional marginal of i^e decays algebraically with a 
power that is a nonincreasing function of b when b is small. It identifies the 
asymptotic behavior of this power as ~ Cst/6 for 6 | 0. Part (e) shows that 
for large b the curve 61— > ^Xm{b) has slope ^(m — 1). 

By Holder's inequality, m 1-^ —Xm{b) is nondecreasing. The system is 
called intermittent of order n if 

(1-47) -Xn{b) < -^Xn+iib) < -^Xn+2{b) < ■■■. 

n n+ 1 n + 2 

(It is shown in [7], Chapter III, that the first of these inequalities implies all 
the subsequent ones.) Thus, for all n > 2 our system is intermittent of order 
n precisely when b G {bn+i-,bn\ (see also Figure 2 in Section 1.1). 

We conjecture that 62 > ^3 > ^4 > • • • [see open problem (A) in Sec- 
tion 1.2]. A partial result in this direction is the following: 

Corollary 1.7. Suppose that a{-,-) is symmetric- 

(a) (m — l)bm — > 2 uniformly in m as G'-^^(0,0) — > 1. 

(b) b2>b3>--->bm w/ien G(2)(0,0) < (m-l)/(m-2). 
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Proof, (a) Obvious from (1.45). 

(b) This follows from (1.45) and G(™)(0,0) > 1. □ 

Claim (a) follows from (1.20) and (1.45), and corresponds to the limit 
when the random walk becomes more and more transient. This includes 
simple random walk on Z*^ with oo. Thus, in this limit all inequalities 
in (1.44) become strict. Claim (b) follows from (1.45). This includes simple 
random walk on with d>3. 

As we will see in Sections 2-3, the representation for bm in (1-43) comes 
from a link with collision local time of random walks. Indeed, let 



(1.48) 




be the total collision local time (in pairs) of the m independent copies of the 
random walk ^. Then we will show that 

(1.49) bra = sup{6 > : (exp[6r(C("^))]) < oo}. 

1.6.2. Regime (II). The next conjecture says that regime (II) is nonempty 
and may therefore be seen as an extension of Theorem 1.6(d). 

Conjecture 1.8. 6* > 62 when a(-,-) is transient. 

Conjecture 1.8 implies that equilibria with stable law tails occur in our 
system for moderate b [recall (1.29)]. We conjecture that the system locally 
dies out at 6^, [see open problem (B) in Section 1.2]. 

In view of (1.49), we may ask whether it is possible to obtain a variational 
characterization for 6*. To that end, let ^ = {C{t))t>o and S,' = iC'{t))t>o be 
two independent copies of the random walk on / with transition kernel a(-, •) 
and jump rate 1, both starting at 0. Let 

(1-50) T{^,a= l{at)=i'it)}dt 

^ 

be their collision local time. Define 

(1.51) 6„ =sup{6>0:E«'(exp[6T(e,C')]) <oo ^a.s.}, 

where we note that {£;«'(exp[6r(^,^')]) < 00} is a tail event for ^. Since 62 
is given by the same formula as (1.51) but with the average taken over both 
^ and ^' [recall (1.49)], we have 6** > 62- The proof of Conjecture 1.8 may 
be achieved by showing that 

(1.52) 6* > 6** and 6** > 62- 

In Section 4, we prove the first inequality and argue in favor of the second 
inequality. A full proof of the latter is deferred to [5] . 
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1.6.3. Regime (III). The last result shows that in regime (III) the system 
gets extinct very rapidly. 

Theorem 1.9 (Cranston, Mountford and Shiga [12]). In regime (III): 
(a) If fi = 6q with G G (0, oo), then 



(1.53) 



lim -logXo(t) = x*(&) 



exists and is constant a.s. 

(b) There exists ah^^ [6=i,,oo) such that 



(1.54) 

(c) 
(1.55) 

exists in (0, oo) 



= 0, forb<K, 
<0, forb>K. 



T log 6 
hm — : — 



X*{b) + -b 



The limit x*(^) is the quenched Lyapunov exponent. Theorem 1.9 states 
that the speed of extinction is exponentially fast above a critical threshold 
6^,. Trivially, 

(1.56) 6*>6,. 

We conjecture that equality holds [see open problem (C) in Section 1.2]. See 
also Figures 1 and 2 in Section 1.1. 



1.6.4. Separation between regimes (II) and (III). The key tool in the 
identification of 6=,, is the notion of Palm distribution of our process X at 
time t. This is the law of the process seen from a "randomly chosen mass" 
drawn at time t. This concept was introduced by Kallenberg [36] in the study 
of branching particle systems with migration. There the idea is to take a 
large box at time t, pick a particle at random from this box (the "tagged 
particle"), shift the origin to the location of this particle, consider the law 
of the shifted configuration, and let the box tend to infinity. Under suitable 
conditions, a limiting law is obtained, which is called the Palm distribution. 
Similarly, in our system the Palm distribution is a size biasing of the original 
distribution according to the "mass" at the origin. The criterion for survival 
versus extinction of the original distribution translates into tightness versus 
divergence of the Palm distribution. 

This criterion is useful for two reasons. First, the size biasing is an easy 
operation. Second, often it is possible to obtain a representation formula 
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for the Palm distribution in terms of a nice Markov process. For instance, 
for branching particle systems the Palm distribution is obtained from an 
independent superposition of the original distribution and a realization of 
the so-called Palm canonical distribution. The latter can be identified as a 
branching random walk with immigration of particles at rate 1 along the 
path of the tagged particle. Fortunately, we can give an explicit represen- 
tation of the Palm distribution of our process X as well, namely, as the 
solution of a system of biased stochastic differential equations (see Section 4 
for details). It turns out that the latter again has a (Feynman-Kac type) 
representation formula for the single components as an expectation over an 
exponential functional of the Brownian motions, the random walk, and an 
additional tagged random walk, with the expectation running over the two 
random walks. 

We will use the Palm distribution to identify 6^,. We will see that, within 
the interval (62,^*), we can distinguish between a regime where the average 
of the Palm distribution over the Brownian motions (i.e., the Palm distribu- 
tion conditioned on the tagged path) is tight as t — > 00 and a regime where it 
diverges. The separation between these two regimes is 6=^,. Within the inter- 
val [6*=K , 6=K ) , we can separate further by conditioning the Palm distribution 
also on the Brownian increments along the tagged path. However, we will 
not pursue this point further, even though it is of interest for a better insight 
into what controls our system. See Birkner [3, 4] for more background. 

We will see in Section 4.1.2 that (1.50) plays an important role in the de- 
scription of the Palm distribution. Equation (4.19) in Section 4.1.2 identifies 
b-f. However, this formula is much harder than the one for 6=^, in (1.51). It 
would be interesting to know whether there exists a characterization of fe* 
in terms of the collision local time of random walks [see open problem (D) 
in Section 1.2], in the same way as for bm in (1.49) and for 6*=,, in (1.51). 



2. Moment calculations. 



2.1. Definition ofbm, bm and bm- For > and m > 2, let 
bm = sup{6 >0:i^e^6o, E^^^iXoD < 00}, 

(2.1) bm. = supfb > : limsupS([Xo(t)]™ | X(0) = 6) < ool, 

I t-*oo J 

6m = sup(6>0:limsup-log^([Xo(t)]™ | X{0) =6) =ol. 

In these definitions the choice of Q is irrelevant as long as > 0, as is evident 
from Lemma 2.1 below. 

In Section 2.2 we derive a representation formula for the solution of (1.2), 
which is due to Shiga [41] and which plays a key role in the present paper. We 
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also derive a self-duality property, which is needed to obtain convergence to 
equilibrium. In Section 2.3 we express the mth moment of a single component 
of our system, at time t, in terms of the collision local time, up to time t, of 
m independent copies of our random walk. In Section 2.4 we prove that 
exists and that bm = bm- In Section 2.5 we prove some basic properties of 
and K^"^^ defined in (1.36) and (1.39). The results in this section will 
be used in Sections 3-4 to prove Theorems 1.4-1.5, 1.6 and 1.9. 

2.2. Representation formula and self-duality. If our process starts in a 
constant initial configuration, then a nice (Feynman-Kac type) representa- 
tion formula is available. This formula will play a key role throughout the 
paper. 

Lemma 2.1. The process {X{t))t>o starting in X{0) = @ can be repre- 
sented as the following functional of the Brownian motions: 



X,(t) = Ge-(i/2)''*£| exp 



^exp 


Vbj 




Jo 



•^0 jei 



(2.2) 

i G /,t > 0, 

where ^ = (^(t))t>o is the random walk on I with transition kernel a{-, ■) and 
jump rate 1, and the expectation is over ^ conditioned on ^(0) = i (£, and W 
are independent). 

Proof. This lemma appears in [41] without proof. We write out the 
proof here, because it will serve us later on. The symbol 1 denotes the 
identity matrix. Note that G enters into (2.2) only as a front factor. 

Fact 1. For all i and t: 

(2.3) X,{t) = e + Vb f y^at-s{i,3)Xj{s)dWj{s) 

Jo . 

with at = exp[t(a — 1)]. 

Proof. Fix i and t. For < s < t, let Yi{s) = J2j at-s{i,j)Xj{s). [The 
infinite sum is finite due to the fact that, by Theorem 1.1(a), [X{t))t>o lives 
in £i defined by (1.4).] Then 



(2.4) dY^{s) 



Y^{l-a)at-s{i,i)Xj{s) 



ds -\-^at-s{i,j)dXj{s). 



From (1.2) we have 
(2.5) dXUs) = J2ia- ^)ihj)Xj{s)ds + VbX,is)dWi{s), 
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which after substitution into (2.4) and cancellation of two terms gives 

(2.6) dYiis) = VbJ2at-sii,j)Xjis)dWjis). 

3 

Integrate both sides from to t, and note that 1^(0) = Yl,j ot(i, j)Xj(O) = 
and Yiii) =J2jO-o{'''ij)Xj{0) =Xi{t), to get the claim. □ 

Fact 2. For all t: 

rt 

(2.7) exp[\/6Zt(t) - \ht\ = 1 + Vbj exp[VbZtis) - ^bs] dZt{s) 
with Zt{s) = Jq J2j l{€(t-r)=i} dWj{r). 

Proof. Fix t. For z G M and s > 0, let f{z,s) = eV^^-(V2)fe^* ^nd put 
g{s) = f{Zt{s), s). Ito's formula gives 

dg{s) = MZt{s),s) dZtis) + y,,{Zt{s),s){dZt{s)f 

(2.8) 

+ fs{Zt{s),s)ds, 

which after cancellation of two terms [because \ fzz + /s = and {dZt{s))'^ = 
ds\ gives 

(2.9) dg{s)=g{s)VbdZt{s). 

Integrate both sides from to t and use that g{Q) = f{Zt{0),0) = 1, to get 
the claim. □ 

The proof of the representation formula in Lemma 2.1 is now completed as 
follows. Let Xi{t) denote the right-hand side of (2.2). Taking the expectation 
over ^ conditioned on ^(0) = i on both sides of (2.7), we get 

e-ix,(t) = ^ + VbJ^ El i|5(i_,)=,}e-ix,(s) dw,{s)^ 

(2.10) ^ 

= l + Vb f y^at-sihj) @-^Xj{s)dWj{s), 
Jo ^ 

where the first equality uses the Markov property of at time t — s. Thus 
we see that Xi{t) satisfies (2.3). Since Xi{0) = Q = Xi{0) for all i, we may 
therefore conclude that Xi{t) = Xi{t) for all i and t, by the strong uniqueness 
of the solution of our system (1.2) [recall Theorem 1.1(a)]. □ 

In addition to the representation formula in Lemma 2.1, we have another 
nice property: our process is self-dual. Let 

(2.11) a*{i,j)=aU,i), hJ^I, 
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be the reflected transition kernel. Let (1.2*) denote (1.2) with a(-, •) replaced 
by a*{-,-). Abbreviate {x,x*) =J2i£iXiXi- 

Lemma 2.2. Let X = {X{t))t>o be the solution of (1.2) starting from 
any X(0) G £i. Let X* = {X* {t))t>o be the solution of (1.2*) starting from 
any X*{0) G £i such that (l,X*(oj) < oo. Then 

(2.12) ^X(g-(xw,x*(o)>) ^ ^x*(g-(x(o),x*w>) > ^ 
Proof. See Cox, Klenke and Perkins [11]. □ 

2.3. Representation of the mth moment in terms of collision local time. 
Let us abbreviate 

(2.13) W = {{Wi{t)hei)t>o 
and write 

(2.14) E{[Xo{t)r I ^(0) = e) = E'^{[Xo{t)r I x{o) = e) 

to display that (1.2) is driven by W. This subsection contains a moment 
calculation in which we use the representation formula of Lemma 2.1 to 
express the right-hand side of (2.14) as the expectation of the exponential 
of b times the collision local time of m independent copies of the random 
walk with transition kernel a{-, ■) and jump rate 1, all starting at 0. 

We begin by checking that the evolution is mean-preserving. This prop- 
erty is evident from (2.3), but its proof will serve as a preparation for the 
calculation of the higher moments. 

Lemma 2.3. E^{Xo{t) \ X{0) = &) = & for all t > 0. 

Proof. Taking the expectation over W in (2.2) and using Fubini's the- 
orem, we have 

E'^{Xo{t)\x{o) = e) 

^^■^^^ 

= ee-(i/2)b*i^U^H/Lxp Vb fY.hat-s)=^}dw,{s 



Since the Brownian motions W are i.i.d. and have independent increments, 
it follows that for any ^: 



2.16) f J2 hi{t-s)=^} dW.i^) ^ W'it), 
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where W' = {W {t))t>o is a single Brownian motion and = denotes equality 
in distribution. Combining (2.15) and (2.16) we arrive at (the expectation 
over ^ being irrelevant) 

(2.17) E^{Xo{t) I X(0) = e) = ee-^^/^^'"E^'{eM'/bW'{t)]). 

Now use that, by Ito's formula, exTp[VbW' (t) — ^bt] is a martingale, to get 
that the r.h.s. of (2.17) equals G. □ 

A version of the above argument will produce the following expression for 
the moments of order m>2. 

Lemma 2.4. Let = {^i, ■ ■ ■ ,Cm) be m independent copies of the ran- 
dom walk with transition kernel a{-,-) and jump rate 1, all starting at 0. 
Then 

(2.18) E^{[Xo{t)r I ^(0) =6) = e"*^«""'(exp[6r(™)(t)]), 
where 

l<k<l<m 

(2.19) ^ 

Tkiit)= dsl{5^(,)=g,(,)}ds, 

is the collision local time {in pairs) up to time t. 

Proof. Similarly as in (2.15) we may use (2.2) to write 

E'^{[Xo{t)r\X{0) = Q) 
(2.20) 

/ / r /-f 



exp 



^/ T.T.h^.it-s)=^}dW,{s) 
■^0 7 — 1 r 



fc=i iei 



Next, let W^'"^^ = {W^, . . . , W^) be m independent Brownian motions. Then 
the analogue of (2.16) reads 

L J2J2h^k{t-s)=i}dWi{s) 

JO I 1 «^ T 



k=l iei 
(2.21) 

rt 



■^0 p=l-ip) q=l 

Here lj{p) [t — s) denotes the indicator of the event that at time t — s the 
components of ^('"^ = {^i, ■ ■ ■ ,im) coincide in p subgroups of sizes j^'P^ = 
(ji) • • • ) Jp) with Ji + • • • + jp = m. The equality in (2.21) again follows from 
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the fact that the Brownian motions W are i.i.d. and have independent in- 
crements. The point to note here is that ah jq random walks in the gth 
coincidence group pick up the same increment of the Brownian motion in 
W at the site where they coincide at time s, and this increment has the 
same distribution as dWg{s). Next, define 

(2.22) T.^,^{t)= j\jM{t-s)ds. 
Then clearly we have 

(2.23) 1^,,) {t - s) dW!j{s) ^ VF^(r.(rt (t)). 
Now combine (2.20)-(2.23) to get 



E'^{[Xo{t)r\X{0) = Q) 



(2.24) 



0mg-(m/2)W£;5(™) / ^ 



exp 



p=l j(p) q=l 



p=i j(p) 



Finally, absorb the term —^ht into the sum by writing \jq{jq — 1) instead 

of \jl [use that Xl^i Ej{p) ^jCp) (*) = * and Eg=i Jg = ^n]. The resulting ex- 
ponent is the same as h times the collision local time in (2.19). □ 

2.4. Convergence to equilibrium and bm = bm- The following important 
facts will be needed later on and will be derived via the representation 
formula in Lemma 2.1 and the self-duality in Lemma 2.2. 

Proposition 2.5. For all a{-, •) transient and all b > 0: 

(a) 1/0 = w-limf_+oo >C(X(t) | X{0) = 0) exists, is shift-invariant and as- 
sociated for all Q G [0, oo). 

(b) is mixing for all Q G [0, oo). 

(c) bm = bm for all m>2. 

(d) limt^oo^([^o(i)]" I ^(0) = e) = ^^e([Xo]'") for all 9 G [0,oo) and 
all m>2. 

Proof, (a) The proof of existence uses Lemma 2.2. If X{0) = Q and 
X*{0) = f, then (2.12) reads 

i?^(e-w*)'^M^(o) = e) 

(2.25) 

= ^^*(e-®^-^W |M(0) = (1,/)) 
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with 

(2.26) M{t) = {l,X*{t)). 

Since {M{t))t>o is a nonnegative martingale [as is obvious from (2.3) with 
the reflected transition kernel], we have that lim^^oo -^^(i) = M{oo) < oo 
W-a.s. for every / [use that ^''^*(M(oo)) < M(0) = (1,/) < oo]. Hence we 
conclude that X{t) converges in law to a limit, which we call vq, given by 

e "<^'-^>i^e(dx) = S^*(e-®*^°° | M(0) = {!,/)) 

(2.27) 

for all / such that (1, /) < oo. 

Because 5q is shift-invariant, so is uq. The fact that is associated follows 
from Cox and Greven [10]. There it is shown that for systems of the type 
in (1.2) — even with a general diffusion term — the evolution preserves the 
associatedness. Since 6e_ is associated, the system is associated at time zero 
and hence at all later times, and the equilibrium inherits this property. 

(b) For < 6 < 62 the mixing property of the equilibrium z^e was proved 
in Cox and Greven [10] via a covariance argument. However, for 62 <i 6 < ^* 
covariances are infinite, and so we must follow a different route. 

The proof uses the exponential duality in Lemma 2.2. We will prove that, 
for all f,g^£i [recall (1.4)] with finite support, 

(2.28) lim ^'^e(e-(^J>e-(^''^fef))=^''e(e-(^,/>)^i'e(e-(^,9>)^ 

where a^g = 9° CTk with the fe-shift acting on /. This implies the mixing 
property, because the Laplace functional determines the distribution. 

Step 1. In order to prove (2.28), we use the self-duality of our process 
and the fact that uq = w-limt^oo <5e5'(t), as follows. Denote by 

(2.29) X*''^ = {X^^\t))t>, = ({X;''^(t)},,^.),>o 

our process with reflected transition kernel o*(-,-) [recall (2.11)] starting 
from initial conflguration h£ £1 with flnite support. Then 

(2.30) = lim E(e-<^W'^+9) | X(0) =G) 

= lim i?(e-<0'^*''^''W)). 

t— >oo 

Observe that, by the linearity of the system, we may use the same Brownian 
motions for X*'^ and X*'^ , which gives us in addition 

(2.31) x*'^+s Ax*J + x*^9_ 
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Hence, in order to verify (2.28), we must investigate the quantity 
(2.32) iimi^(e-(e,x*./(t)+x*-.^(i)>) 



and show that it factorizes in the hmit as ||A;|| — > oo. 
Next, note that 

(2.33) m,X*'fmt>o and ((6, ^*'"'=^(t)))t>o 

are (continuous-path square-integrable) nonnegative martingales. In partic- 
ular, their limit as t ^ oo exists by the martingale convergence theorem. 
Their covariation over the time interval [0, oo) is given by 

C{f,akg)= / ds{X*'fis),X*'^'^^{s)) 

In 

(2.34) 



/■oo 

/ dsY^X:'f{s)Xr^^{s). 
Jo j^l 



Due to these structural properties, we know that if 

(2.35) C{f,ai^g) in probability as — > oo, 

then the two martingales in (2.33) become independent as ||A;|| — > oo. Con- 
sequently, the two random variables 

(2.36) lim(G,X*'^(t)) and lim (G, X*'^'=»(t)) 

t—^00 t—^oa 

also become independent as ||A;|| — > oo, which proves (2.28) via (2.30) and 
(2.31). In order to prove (2.35), observe that, by the linearity of the system, 
it suffices to verify (2.28) for the special case where / and g are indicators 
of a single site in /, say, p and q, respectively. 

Step 2. Let ^ and ^' be two independent random walks with transition 
kernel a*(-, •) and jump rate 1, both starting in i G /. Then, for / = Ijp} and 
g = l^gy, it follows from Lemma 2.1 that 

x*'f{s)xr''\s) 

e-'^Elf(^ms)=p,eis)=q + k} 

rs 

Vb du [l{^(-s — u)= m} 
^0 me/ 

+ l{i'{s-u)=m}] dW^{u) 



(2.37) 

X exp 



BRANCHING PROCESS IN A RANDOM ENVIRONMENT 



25 



Reversing time, we may start ^ in p and ^' in q + and give them 
transition kernel a(-, •) and jump rate 1. Then 

Vb Tdn V[l{^(n) = m} 

Jo 



(2.38) 



X exp 



mei 



+ l{C'{u)=m}]dWm{u) 



and so 



(2.39) 



X exp 



+ l{e'(u)=m}] dWm{u) 



Next, for any s we have 



(2.40) 



where 



(2.41) 



{X*^f{s),X*''^^3{s))=Y^Xl^i{s)Xr^\s) 
= iMi(s) + iM2(s), 



^exp 


2\/6 / 




Jo 



^exp 


2Vb [ 




Jo 



2Vb f duY l{C{u) = m} dWm{u) 
^0 me/ 
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Below we will show that 

(2.42) / dsMi{s) <oo W^-a.s. for i = 1, 2. 

Jo 

Assuming (2.42), we pick T > and estimate, with the help of (2.40), 



(2-43) 

< / ds{X*'f{s),X*'''^3{s)) + ^ r [Mi{s) + M2{s)]. 



T 

By (2.42) and the fact that the law of (M2(s))s>o is independent of k, it 
now suffices to show that 

/^(is(X*'^(s),X*'"'=5(s)) ^0 
(2.44) ^° 

in probability as ^ oo for any T > 0. 

Step 3. To prove (2.44), we return to (2.39). Write 

miu) =m} + l{e'(n) = m}] dWmiu) 

= l{ay')=^'{u)=m}2dWm{u) 

(2.45) 

+ Y ^'{u)} dWm{u) 

+ E l{e(n)7^m = e'(n)}dT^^(n). 

Since the three terms in the right-hand side of (2.45) involve disjoint time 
intervals and the Wm's have independent increments, it follows from (2.45) 
that, for any ^, 

Vb r duY\ =m} + l{i'{u) = m}] dWm{u) 

(2.46) 

^ Wi{AhTs{W)) + W2{h[s - r,(e,0]) + - r.(e,0]), 

where Wi,W2,W3 are three independent Brownian motions, and Ts{^,^') = 
Jq dul{S,{u) = £,'{u)} is the collision local time of ^ and ^' up to time s. 
By combining (2.39) and (2.46), taking the expectation over W (i.e., over 
VFi,VF2,VF3) and using Fubini's theorem, we get 

rT 
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(2.47) =^'^ii;J;«;^(l{e(.) =e'(.)}e^^»(«'«')) 

Clearly, for fixed T the right-hand side tends to zero as ^ oo, because 
Tt{C,£,') 1^ T and Tt{C,(,') ^ in probability with respect to as \\k\\ 
oo for any fixed T. 

Step 4. It remains to prove (2.42), which goes as follows. Let 

(2.48) (M(s)),>o withM(s) = (l,X*'^(s)). 

This is a (continuous-path square-integrable) nonnegative martingale start- 
ing from a strictly positive and finite value (because / = Ijp} has finite 
support). From the dual of (1.2) (recall Lemma 2.2), we have 

dM{s) = J2dX:'^{s) 

(2.49) =Y.a*{i,j)[xf{s)-X:'f{s)]ds 

+Y,^h[x:^f{s)Ydw.{s). 

The first term in the right-hand side is zero because a* (•, •) is doubly stochas- 
tic (being a random walk transition kernel [recall (1.9)] and {1, X*'^ (s)) < oo. 
Hence 

(2.50) M{s) = W{t{s)) 
with 

(2.51) r(.)= f duhY^lXl'f {u)f 

and W a Brownian motion adapted to the filtration of X*'f . By the mar- 
tingale convergence theorem, we have 

(2.52) lim Mis) = M(oo) < oo, W-a.s. 

S — s-OO 

[use that E^*'\m{oo)) < M(0) = (1,/) < oo]. Combining this with (2.50), 
we conclude that 

(2.53) lim t(s) = t(oo) < oo, W-a.s. 

S — >CXD 

This completes the proof of (2.42), hence of (2.35), and therefore also of 
(2.28). □ 
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(c) Fatou's lemma in combination with part (a) shows that 

(2.54) limmfE{[Xo{t)r \ ^(0) = 6) > ^'^^ ([^o]™)- 

Hence bm < bm - The converse is proved as follows. Assume that z^e 7^ ^o- 
Define the jn-point correlation function in equilibrium, 

(2.55) /(ii,...,j^)=£^«(^nx,-^^, ju...,jmel, 

where the indices need not be distinct. Since i^q is associated and shift- 
invariant [which was proved in part (a)], we have 

(2.56) < fin, ...,jm)< /(o, . . . , 0) = E'^^iiXoD. 

Moreover, from the equilibrium property of i^e we deduce that, for any t>0, 
E""" (ft ^^•.) = / Mdx^E'^iiXoit)]"^ I X{0) =x) 

(2.57) 

= i?«^™^(exp[6T„(t)]/(ei(t),...,r (*))), 

where the last line follows after substituting the representation formula in 
Lemma 2.1 [with an arbitrary initial condition ^(0)] and doing a calculation 
similar to the one in the proof of Lemma 2.4. Passing to the limit t ^ oo in 

(2.57) , we get, with the help of (2.56), that 

(2.58) ^'^^([Xo]'") <oo =^ (exp[6Tm(oo)]) < oo. 
Hence bm>bm- 

(d) We need to show that (2.54) is an equality. This is trivial when the 
right-hand side of (2.54) is infinite. Therefore, assume that uq ^ 6o and 
E'^^ {[Xoy^) < oo. By applying the mixing property of z^e [which was proved 
in part (b)] to (2.57), we get 

(2.59) E'^^iiXoD = e^E^'"'\eMbTmioo)]), 

where we use (2.58), dominated convergence and the fact that ,^(*")-a.s. all 
m random walks move apart as t — > oo by transience. Moreover, by passing 
to the limit t ^ oo in (2.18), we have 

(2.60) lim E{[Xo{t)r \ X{0) = 9) = 6™ S^'"' (exp[6r,„(oo)]). 
Combine (2.59)-(2.60) to get the claim. 
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2.5. Properties of G^™^ and K^™'\ This section lists a number of ele- 
mentary facts, many of which use basic random walk theory as explained in 
Spitzer [43]. 

Recall (1.36) and (1.39). We begin with the following statement: 

Lemma 2.6. Ifa{-,-) is symmetric and transient, then G^"^\-, ■) is strongly 
transient for all m>3, that is, 

(2.61) sup G^"'\x,z)G^"'\z,y)< oo, m>3. 

Proof. Let 

P/"*) {x,y)= P^'"^' (7?(™) {t) = y\ r/(-) (0) = x) , y E /(™) , t > 0, 

(2.62) 

be the transition probabilities of the differences random walk 77*- '"^ = (r^^™) (t))t>o 
defined in (1.33). We have 

POO , , 

(2.63) G^"'\x,y) = m dtPr\x,y), x,y£l^"'l 

Jo 

Compute 

G^^\x,z)G'^^Hz,y) 

) /"OO 

ds^ \ ds2 Y PtH^,^)Pt^iz,y) 

PCSD /'OO 

(2.64) =/ dsj ds2PtLix,y) 

Jo Jo 



roo 



dt tPt^'^\x,y) 

< r dt tpj:"'\o,o). 

Jo 



Note that 

(2.65) P,^"'\0,0) = Y[Pti0,i)r 

with Pt{i,j), i,j(zl, t> 0, the transition probabilities of a single random 
walk. Because the single random walk is symmetric and has exponential 
jump times (with mean 1), we have 

(i) Pt{i,j)<Pt{o,o)yi,jGi,t>o, 

(2.66) 

(ii) 1 1-^ Pt(0, 0) is nonincreasing, 
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as is easily seen from the Fourier representation of Pt{i,j); that is, 



(2.67) 

= {Iny^ I dX cos(i - j, A)e-*[i"^(^)l 

with A{X) = X^ig/ a(0, i) cos(i, A) < 1, AG [—tt,ttY, and (•,•) the inner prod- 
uct on E'='. Via (2.66) (i), (2.65) gives 

(2.68) Pt\0,0) < [Pi(0,0)]™-2^Pi(0,i)Pi(i,0) = [Pt{0,0)r-'P2t{0,0), 
which via (2.65) (ii) yields 

/"oo /"oo i''2it 

/ dttP^'^\<d,Q)<\ dt[Pt{0,0)r-^ dsPs{0,0) 
Jo Jo Jo 



(2.69) 



dt [Pt{0,0)r~^)[l"° dsPMO)^ 



The right-hand side is finite by the transience of the single random walk. 
□ 

We next look at K^™) defined in (1.39). 

Proposition 2.7. Suppose that o(-, •) is symmetric and transient. Then, 
for all m > 2, K^^^ is a self-adjoint, positive and bounded operator on 

Proof. The symmetry of follows from the symmetry of G^"^\ 

Since K^"^^ is defined everywhere on i'^{S^"^'^), it therefore is self-adjoint. 
The Fourier representation of G^^^ reads 

oo 

G(™)(x,y) = ^[a(-)r(x,y) 

71=0 

(2-70) 

= dA e^(^-^'^)[l- A('^)(A)]"\ x,y£l^"'\ 

with J^'") = ([-TT, 7r)^)(i/2)m(m-i) ^ (.^ .) ^j^g ^^^^^ product on (Mrf)(i/2)m(™-i) ^ 
and A^"''^{X)=J2xei(m) a^™) (0, x) cos(x, A) < 1, A G /(""^ It follows that 

(2.71) 

dA [1-A(™)(A)]-^ V e^(^'^V(2;)\/tt(™)(^) 
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with (•, •) the inner product on i'^{S^"^^). This proves the positivity of K^"^\ 
To prove the boundedness of K^^\ we consider the relation 



(2.72) 



livx^l 



: sup {^l,K^'^^^i) 

{m>m>=i 



with II • II denoting the operator norm on £^(S'("^^). Apply Cauchy-Schwarz 
twice, to obtain 



x,yeS("^') 

,ygS(™) 



< E 



(2.73) 



< 



J2 /i'(x)G(-)(x,y)tt('")(y) 

1/2 

^ «M(^)G'M(x,y)/.2(y) 



1/2 



^ /i2(x)G(-)(x,y)ttM(y) 



.1/2 



< 



sup ^ G(-)(x,y)S(™)(y), 



where in the last line we use the symmetry of G^™-). The last sum is equal 
to the average total collision local time (in pairs) of the m walks when their 
differences start in x. Clearly, the supremum is taken at x = 0, and equals 
ljM(o)G(2)(0,0), because G(2)(0,0) is the average collision local time for 
each pair. Hence 



(2.74) II^KMf < tt("^)(0)G(2)(0,0) < oo. 

But, by the self-adjointness and positivity of K^'^\ we have (see [39], Chap- 
ter 12) 



(2.75) II Vi^M f = II ^ spec(i^('")) VnGN 

with spec(-) denoting the spectral radius in £'^{S^"^^). □ 
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3. Variational representations. In Section 3.1 we identify bm in terms of 
a variational formula. In Section 3.2 we prove that the mth moment diverges 
at 6 = 6m. In Section 3.3 we calculate the exponential growth rate of the mth 
moment and prove that bm = bm- In Section 3.4 we study the m-dependence 
of bra- In Section 3.5 we collect the results and prove Theorem 1.6. 



3.1. Variational formula for h, 



Proposition 3.1. Suppose that a(-, •) is symmetric and transient. Then 
bm = rn/Xm with 



Proof. The proof comes in several steps. Throughout the proof we 
assume that 6jj('")(0) < m. 



Step 1. Recall the definition of bm in (2.1) as well as the identity in (2.18). 
We begin by deriving a criterion for the property E{exp[bT^"^\oo)]) < oo in 
terms of the discrete-time random walk 

(3.2) ryM,0 = (^M,0(,)).^^^ 

embedded in the continuous-time random walk ry^™) = {r]^"^\t))t>o defined 
in (1.33). To that end we perform the expectation over the jump times of 
fji^) ^ which are independent of r/^"^)'®. Indeed, let 

(3.3) (cJi)igNo 

be the successive discrete times at which ry'-'"^'® visits S^"^^ (put ctq = 0), 
and let M be its total number of visits to S'(''") (which is random but finite 
a.s. by transience). Each visit to S^"^^ lasts a time r that is exponentially 
distributed with mean — . Define 

m 

(3.4) n^-^)^\x)=E{eMH^'^Hx)T]) = -—^^^, xe5(-). 

Then we have 

(3.5) £;(exp[6T(™)(oo)]) =i?(^n°^™^'''(^^"^'®(^^))) ■ 
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Step 2. In order to analyze the right-hand side of (3.5), we introduce 
the transition kernel of the Markov chain on S^"^^ obtained by observing 
fji"^)'® only when it visits S^"^\ which we denote by P^"^\-, •). Since a(-, •) is 
symmetric, so is (•,•). By transience, this transition kernel is defective: 

(The first line says that escape from S^"^^ is possible only when all walks are 
disjoint except one pair. This is because only one walk moves at a time.) In 
terms of P^"^\-, •) we can write 

oo / n \ 



n 



X 



Ki=0 



Define the matrix 



(3.8) QM'^(x,y) = y^nM.''(x)p(™)(x,y)y^n(™).''(?/), x,yeS^"'\ 

With this notation we can write, combining (3.5) and (3.7)-(3.8), 

oo 

(3.9) E(exp[6T(™)(oo)]) = C™,ft^(5o,[Q('")'T^m,) 

n=0 

with 



and 



Cm,b = m/y'(m-6SM(0))(m-6) 
ii^(-) = [l-P("^)(-,5(-))]%.)(.)=i|. 



The front factor, which arises from the endpoint in the second sum in (3.7), 
is harmless. 

Step 3. Note the following: 

Lemma 3.2. Q(™)'''(-,-) is an irreducible, aperiodic, nonnegative and 
symmetric matrix. As an operator acting on £^(5^™)) it is self-adjoint and 
bounded. 
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Proof. Because is bounded, •) inherits these properties 

from P^"^\-,-). The irreducibihty of P^"^\-,-) is inherited from the irre- 
ducibihty of a(-,-) assumed in (1.9). The aperiodicity of P^'^\-,-) follows 
from the fact that P^"*) {x, x)>0 for some x G S^""^ with tt^"*) {x) = l. □ 

Define 

(3.10) xn.{b) = lim ilog[Q(™)'T(x,y), x,yeS^^\ 

n — *oo 72 

Under the properties stated in Lemma 3.2, this hmit exists, is in M and is 
the same for ah x,y G S^™) (see Vere- Jones [45]). Moreover, 



^[Q(-)'T(0,0)<oo 



(3.11) 



n=0 



£[Q(™)'T(a;,y) <oo yx,y£S^"'\ 



n=Q 



This leads to 



(3.12) 



Xm{b)>0 =^ X][Q^'"^'T(0,0) = oo, 



n=0 

oo 



Xrn{b)<0 =^ ^[Q(-)'T(0,0)<oo. 



n=0 

From (2.1), (2.18), (3.5), (3.9) and (3.12) we see that bm is the solution 
of the equation Xmib) = 0. At the end of Section 3.2 the case b = bm will be 
included in the top line of (3.12). 

Step 4. Next, at b = bm we have the following: 

Lemma 3.3. 

(3.13) sup 7 = 1. 

Proof. Since is a bounded operator on ^^(S^"^)), the function 

(i/,Q('")'S-i/) is continuous on ^^(sM). Since ^i(5(™)) C ^^(^M) is 
dense, it suffices to prove that 

(3.14) sup (z^,Q("^)'^™i/) = 1. 

j,g£2(5{m)) 

{u,u)=l 
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[< 1]: First consider v with finite support. Suppose that (z^, Q^™)'^™!/) > 
1 + e for some e > 0. Then, by the spectral theorem and Jensen's inequality, 



(3.15) 



with ii^jy^i, the spectral measure associated with v. Clearly this contradicts (3.10) with 
Xm(bm) = 0, and so {y, Q^"^^'^"^!/) < 1 for u with finite support. Since the z/'s 
with finite support are dense in it follows that the supremum is 

< 1. 

[> 1]: Suppose that the supremum is < 1 — e for some e > 0. Then the 
spectrum of is contained in {—oo, 1 — e]. Estimate 

o< (5o,[g("^)'"'"f"+'5o)= r\^''-^UEs,M{\) 



(3.16) 

"1— e 



oo 



</ A2«+idi?5„,5o(A)<(l-e)'"+^ VnEN. 







But again this contradicts (3.10) with Xm{bm) = 0. Hence the supremum is 
> 1. □ 



Step 5. Putting ^ = vTlK™7^^z^, we may rewrite (3.13) as [recall (3.4)] 



(3.17) 



Therefore 



1 = sup , , , 

(/i,/i)-(M,(l-P("-))/i) ^ 



(3.18) hra = ^ withA„= sup ^/^'S^"^/^) 



Am MS^HSf-)) (1 - ^'('"))/^> ' 

where the denominator is strictly positive because P^"*) is irreducible. Let 

oo 

(3.19) G('^) = ^[pM]" = (i-pM)-\ 

n=0 

Because P^'") has spectral radius < 1 (due to the fact that 5^™) is a uni- 
formly transient set), we know that 1 — P^^' is one-to-one on ^i(5('")) and 
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that G(™) is a bounded operator on £^{S^"^''). Therefore we can transform 
(3.18) via the change of variables fj, = G^"^'^ p: 

(3-20) Xm= sup , . , • 



Finally, putting p= -\/tj("^)C and using that G^"^\x,y) = G^'^\x,y) for all 
x,y £ 5^"*^ by the definition of P^™-), we get the formula in Proposition 3.1. 
□ 



3.2. The mth moment atb = hm- The case h = bm can be included in the 
top line of (3.12) when 1 is the largest £^-eigenvalue of Q'-™^''''". Therefore 
we next consider the eigenvalue equation 

(3.21) z.qM.5™=z., i/E^i(5(™)), i/>0. 

Lemma 3.4. Suppose that a(-,-) is symmetric and transient. Ifbm-i > 
bm, then (3.21) has a solution. 

Proof. The idea is to use the notion of a quasi-stationary distribution. 



Step 1. Consider the matrix 
(3.22) 

This is an irreducible defective probability kernel on S^'^\ By introducing 
a cemetery state 9, we can extend q(™)'^»"'® to a nondefective probability 
kernel on S^'") U {5}. Let (Zn)neNo 

denote the corresponding Markov chain 

starting in 0, and let 

(3.23) Vn = C{Zn\Zn^d), n G Nq. 

If we manage to show that {V denotes the set of probability measures) 

(3.24) limz.„ = i/oo inP(5(™)), 

n—*oo 

then, because 



(3.25) Vn+l 
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we get that 

(3.26) i/ooQ^™^''^'"'® = Aooi^oo, Aoo > 0, z^oo > 0, 

with Aoo = (I'ooQ^™^''''"'®, 1) = the probabihty of no defection to d 

(per step) in the quasi-stationary distribution z^oo- Hence Voo solves (3.21). 

Step 2. To prove (3.24), we use a criterion in Ferrari, Kesten and Mar- 
tinez [23], Theorem 1, according to which it is enough to prove that there 
exist (5 > and D < oo (depending on m) such that 

(3.27) Po{To>n\Ta>n)<De~^'' Vn G No 

with Pq the law of {Zn)n£No given Zq = 0, and tq,tq the first hitting times 
of 0,9 (time zero excluded). 

For K C{l,...,m} with < | A'| < m, let 

Vk = {x G S^"^^ : for site x there exist i ^ K,j £ 

(3.28) 

such that walks i and j coincide} 

with = {1, . . . \ K (recall from Section 1.6.1 that each x S /'-'"^ cor- 
responds to a certain intersection order of the m random walks). We will 
prove that there exist 5 > and Dk < oo such that 

(3.29) Po{TVj,>n\Td>n)<DKe-^'' VnGNo, 

where ry^ is the first hitting time of Vk (time zero excluded). Since 

(3.30) fl Vk = {0}, 

0<\K\<m 

we have {tq > n} C Uo<|i^|<m{''"Vif > "'^l) a-i^d so (3.29) implies (3.27) with 
D = J:kDk. 

To prove (3.29), write 

Po{rvK>'^\'Td>'n) = N/D with 

N= Q("^)'^™'^(0,yi)Q('")''^'-®(yi,y2) x ••• 

j/i,...,j/„es('")\yx 

(3.31) xQ(-)'"^'"'«(y„_i,y„), 

D= Yl Q^'"^''"^'®(0,yi)g^'")'''"'^(yi,y2) x ••• 
xQ('")'"''"'®(y„-i,y„) 
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We may drop the ® because the normahzation factor in (3.22) cancels out. 
After that the denominator in (3.31) equals 

{So, > {do, [Q('")'"^"T'5o) = [XM + 0(1)]" = exp[o(n)] 

(3.32) 

because Xm{bm) = 0. It therefore suffices to prove that the numerator in 
(3.31) satisfies the exponential bound in (3.29). 
Now, on S'('") \ Vk we have 

(3.33) r('")(a„)=r('")'^(a„)+r('")'^'(a„) V?iGN 

with Un the time of the nth visit to [recah (3.3)] and r('")'^(CTn) the 

total collision local time (in pairs) up to time cj„ of the walks indexed by K, 
and similarly for K'^. Therefore, retracing the calculations in Steps 1 and 2 
of the proof of Proposition 3.1, we find that 

numerator (3.31) < i?o(exp[6^{r('")'^(aO + r("^)'^^(a„)}]l|,„<oo}). 

(3.34) 

The inequality comes from using (3.33) and afterward dropping the restric- 
tion to S^"^) \ Vk- Next, apply Holder's inequality to get, for e > 0, 

numerator (3.31) 

< i?o(exp[5^{rM'^(cx)) +r(-)'^^(oo)}]l|^„<„o}) 

(3.35) 

< i?o(exp[(l + e)C{T(™)'^(oo) + T^^)^^\<x,)]]f^^^'^ 
xPo(an<oo)^/('+^). 

The expectation in the right-hand side factors because r("^)'^(oo) and T(™)'^''(oo) 
are independent, and each factor is finite when e is picked so small that 
(1 + e)bm < bm-i, because \K\, \K^\ <m — l. On the other hand, the proba- 
bility in the right-hand side tends to zero exponentially fast with n because 
5"^™) is a uniformly transient set. □ 

By (2.1), we have E^(^{[Xo]'^) < oo for 6 < 6„ and ^^©([Xq]'") = oo for 
b > bm- With the help of Lemma 3.4 we can now include b = bm- 

Lemma 3.5. E""® {[Xo]^) = oo at b = bm. 

Proof. By Proposition 2.5(b), we have bm = bm- Suppose first that 
bm~i > bm- Then, by Lemma 3.4, 

(3.36) z.Q^'")'^'" = 1/ u££^ (5(™) ) , for some > 0, 
and hence 

oo oo 

(3.37) ^ ^ u{x)[Q^^'^M^{x,y) = Y,u{y) = ^ ^yeS^^l 

n=0^.g5(™) n=0 
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By the irreducibility of Q^™)'''™, this imphes that 

oo 

(3.38) ^[Q('")''''"]"(0,0) = oo, 

n=0 

which shows, via (3.5) and (3.9), that 

(3.39) E«*'"'(exp[6^r(™)(oo)]) = cx). 

Now use Proposition 2.5(d) to obtain from (3.39) that E"'-^ {[XqY'') = cx) at 

b = fern- 
Suppose next that ftm—i = bm, but bm'-i > bm for some m' < m. We may 

assume that m' is the largest such index. Then bm'-i > bm', and so the above 

argument tells us that E'^^{[Xq]'^ ) = oo at 6 = 6^'- But because bm' = bm 

and m' < m, we again get E"^ {[Xq]"^) = od at 6 = bm- 

Finally, if bm' = bm for all m' < m, then 62 = ^m- However, at 6 = 62 we 

have -E''®([Xo]^) =00, as is easily seen from Lemma 2.4, Proposition 2.5(d) 

and (3.9), because 5^2) = {0}. Therefore once again E'^^{[Xo]'^) = 00 at 

b = bm- □ 



3.3. Growth rate of the mth moment and bm = bm- In this section we 
show that 

(3.40) hm ^logii;«""'(exp[6r(™)(t)]) = x™(6) 

exists and can be expressed in terms of a variational problem. We will analyze 
this variational problem and show that 

(3.41) 6^ = sup{6>0:xm(6) = 0} 

[recall (2.1) and (2.18)] coincides with bm- Together with Proposition 2.5(c) 
this will show that all three critical values in (2.1) coincide. 

In order to pose the problem in a form suitable for a large deviation 
analysis, we recall the definition of the differences random walk {ri^^\t))t>o 
in (1.33) and the collision function jj^™^ in (1.37). Using (2.19), we have the 
identity 

(3.42) r('")(t)= /V™'^(r/^"'^(s))(is. 

Jo 

Next, we define the empirical measure 

(3.43) ^*"^ = ii* 
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Lemma 3.6. Suppose that a(-,-) is symmetric. Then (L|™^)f>o satisfies 
the weak large deviation principle on with rate function 

(3.44) j('")(zy) = (z/i/2,m(l-o(™))i/^/2^, veV{I^'^^), 

where a^"^\-,-) is the transition kernel defined in (1.34). is bounded on 
£2(/M) and when restricted to (/(")) is continuous in the -topology. 

Proof. See Deuschel and Stroock [21], Section 4.2. The rate function is 
given by (3.44) because m{a^"^^ — 1) is the generator of the Markov process 
{r]^"^\t))t>o and a^"^\-, •) is symmetric [recall (1.32)]. The latter is crucial for 
having the explicit formula in (3.44). The boundedness of J^"^^ is obvious. 
The continuity of J^"^^ follows from the fact that m(l — a^™-*) is a bounded 
operator from ^^(/(™)) into in combination with the fact that if 

f„ — > as n ^ oo in -f-'^-norm, then z/.y^ u^^^ in ^^-norm. To see the latter, 
write 



(3.45) - v^lX2 = E(V^-(^) - V^(^) ) = 2 -'i^Y.^^n[x)v{x). 

X X 

liun^ I' in ^^-norm, then Fatou's lemma gives liminf^^oo X^x V i^n{x)v{x) > 
J2x ^{x) = 1- Hence lim„^oo 11^'^^ — ^^^"^Wi = 0- ^ 

Lemma 3.6 leads us to the following identification: 
Lemma 3.7. 

(3.46) xm{b)= sup jM(z.)}. 

Proof. For ease of notation we drop the superscript (m). We cannot 
apply Varadhan's lemma directly to (3.40)-(3.43), since we only have a 
weak large deviation principle. This problem can be handled via a standard 
compactification argument, as follows. 

Let (??^(i))i>0) ('?Ar(0)<>o be the differences random walks obtained by 
wrapping (r?(t))t>o around the torus Aat = ([-iV, AT)*^ n /)(^/2)'"(™-^) (recah 
that / = Z'^), respectively, by killing it when it hits dA]\!, the boundary of 
Ajv. Let Tjy(i),T^(t) be the quantities corresponding to T{t) in (3.42) for 
these two processes. Then, by the large deviation principle in Lemma 3.6 
restricted to A^r, we have for every N: 

lim ^logE(exp[5r+(t)]) = 5+, 

(3.47) 

hm ilog^(exp[5r^(t)]) = 5]^, 

t — >oo t 
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with 

Sjf= sup {b{i^N,\!,) - Jn{i^n)}, 

(3.48) 

5]^= sup {bivN^fl-JNil^N)}- 
i^jv{9Ajv)=0 

Here, Jn is the analogue of J in (3.44) restricted to A^r, that is, 

(3.49) JN{vN) = {vTM^-aN)vT)^ meViAN), 
with aN{-,-) the periodized transition kernel 

(3.50) aN{x,y)= ^ a(x,2/ + 2:), x,yeAN. 

2e{2Af)/(™) 

Now, obviously 

(3.51) T+{t)>Tit)>T^{t), 

because the intersection local time increases by wrapping and decreases by 
killing. Consequently, 

(3.52) S+>S>S^ 

with 5 the right-hand side of (3.46). Hence, to prove (3.40) and (3.46) it 
suffices to show that 

(3.53) liminf5^>5, lim sup 5]y < S*. 

N^oo N^oo 

Lower bound. It suffices to show that for every z/ E 'P(/'-'^^) there exists 
a sequence (i^n), with i^n £ 'P(AAr) and ^^{dAN) = for each N, such that 

(3.54) W-liniAT^oo l^N = l^i hm JAr(z^Ar) = Jii^)- 
Indeed, with the help of Fatou's lemma this gives 

(3.55) liminf 5^ > liminf{6(z^Ar, ft) — JNi^^N)} > b{i^, jj) — J{i^), 

and so we get the first half of (3.53) after taking the supremum over u 
afterward. 

For the given v, the sequence (i^n) is chosen as follows. Put vn{x) = v{x) 
for all X € An \ (OAn U {0}) and i^7v(0) = i^{0) + Z]x^Ajvu9Ajv ^i^) - Then the 
first half of (3.54) is obvious. For the second half, since J is continuous in 
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the £i-topology it suffices to show that lim]\f^oo[JN{i^N) — J{^n)] = 0. To 
that end, we estimate 

< J{i^n) - Jn{vn) 

/ 1/2 , . l/2\ 

(3.56) <m ^ UN{x){aN - a){x,y) 

<m sup ^ {aM - a){x,y) 

<msup ^ a{x,y) = 5N, 
\\y-x\\^>N 

where in the third hne we use the Cauchy-Schwarz inequality and the sym- 
metry of ttiv and a, and in the fifth hne we use (3.50) and the shift-invariance 
of a(-, •) (note that on the sublattice with spacing 2A^ containing y G A^r, the 
site closest to x is captured by the supremum over x G A^v and the sum over 
y £ ji^) with \\y — x\\oo < N). Obviously, liniAr^oo = 0, which completes 
the proof of the first half of (3.53). 

Upper bound. Estimate, with the help of (3.56), 

S^= sup {b{iyN4) - Jn{i^n)} 

i^ivSPCAjv) 

(3.57) < sup {b{i^N4) - Ji^N)} + 
<S + Sn. 

Let — > oo to obtain the second half of (3.53). □ 

It follows from (3.46) that b Xm{b) is nondecreasing and convex on 
[0, oo). Since it is finite, it is also continuous on [0, oo). Obviously, Xm(0) = 0. 
The critical value bm is given by (3.41). It further follows from (3.46) that 

(3.58) bt'Ho) - J^^'HSo) < Xm{b) < 6tt('")(0), 

where J^'^HSq) = m{l - a(™)(0,0)) = m [use (3.44) and (1.34), and recall 
that a(0,0) = as assumed below (1.9)]. 

Proposition 3.8. Suppose that a(-, •) is symmetric and transient. Then 
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Proof. Changing variables in (3.46) by putting v = , we get 

(3.59) xm{b)= sup f(-)(/.) 

with 

(3.60) f('")(^) = tt(™)) - m(^, (1 - a(™))/i). 
Define 

(3.61) Am = sup 



/.e^2(/{-)) (/^^Cl-a^""^)/^)' 

where the denominator is strictly positive because a^*"^(-, •) is irreducible. It 
follows from (3.59)-(3.60) that 

h\m >m =^ Xrn{b) > 0, 

(3.62) 

bXm<m =^ Xmib)=0. 

Hence bm = m/^m- Change variables in (3.61) by putting = G^^^ p. Then, 
because G^"') = (1 - a(™))-\ we get 

(3-63) A„> sup — , ^( ^ , — . 

Here the inequality arises because we restrict the support of p to S^'^'^ : G^"^^ p 
does not run through all of £^(5'-™'^). Note here that, because G*-'"^ is strongly 
transient by Lemma 2.6, p £ ^^5^™)) implies that p = G^"''^ p G £2(j(m))_ 

Putting p = ^fiy^C with C £ t{S^'^^), we obtain 
(3.64) A„> sup ^)f^J)X - 

Comparing with (3.1), we have thus found that Am > Am,- Consequently, 
bm ^bm- However, it is obvious from (2.1) that 6m ^ bm- Hence we get 
bm = bm and Am = Am- □ 

3.4. Analysis of Tfi I — > bm. • The remaining steps in this section are the 
following three propositions. 

Proposition 3.9. Suppose that a(-, •) is symmetric and transient. Then 
bm = bm = bm = Y~ ^^^^ "^"^ Spectral radius of K'^'^^ in £^(5*^™)). 
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Proof. It follows from Propositions 2.5(b), 3.1 and 3.8 that bm = bm = 
bm = with Am given by the right-hand side of (3.1). We saw in the proof 
of Lemma 3.3 that the supremum may be taken over u G £'^{S'^"^^). Putting 
fi = VK(^iy in (3.1), we see that = . Now recall (2.75). □ 

Proposition 3.10. Suppose that a{-, •) is symmetric and transient. Then 
the spectral radius of K^'"^^ is an eigenvalue if and only if bm.^i > bm- 

Proof. The fact that 6m-i > bm implies that the spectral radius of 
is an eigenvalue follows from Lemma 3.4 via the change of variables 
in Step 5 in the proof of Proposition 3.1. Indeed, Lemma 3.4 shows that the 
supremum in (3.1) is then actually attained. 

To get the reverse, we use Lemma 3.4 and Proposition 3.9. The idea is to 
imbed the variational problem with label m — 1 into the one with label m. 
To that end, consider m random walks but ignore the collision local time 
the mth random walk has with the others. By repeating the large deviation 
analysis in Section 3.3, we get a formula for Xm-iib) of the same form as in 
(3.46) but with the first term (i^, tt^™^) replaced by 

(3.65) (i/,tt('")'("'-^)), 
where [cf. with (1.37)] 

l<p<q<m~-l 

(3.66) 

In this way we can represent Xm{b) and Xm--i{b) via a variational problem 
on the same space: 

Xm{b)= sup {6(z.,tt(-))-j(-)(i.)}, 

(3.67) 

X^-i{b)= sup {6(z.,tt(™)-(™-i))- j(-)(z/)}. 

Now consider these two variational problems at b = bm- Suppose that K^^^ 
has a maximal eigenvalue. Then this eigenvalue is unique, and so is its 
corresponding eigenvector. Consequently, the first variational problem has 
a unique maximizer, say P G V{I^"^^), which is strictly positive. For e G 
(0,z7(0)), put 



(3.68) 



^, = {i.GP(/('")):|z.(0)-P(0)|<e}, 
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and let 

A, = sup J(")(i/)}, 

(3.69) 

B, = sup{6„(i.,S('-)'('"-^)) - J^'^Xu)]. 
Since v is the unique maximizer, we have 

(3.70) ^ = Xm{hm)> sup {hm{v,i'^^)-^'^\v))>A,. 

Since tj^'")(0) > t|("')'("'-^Ho), we have 

(3.71) = XmM>snp{bmWA^"'^)-J^"'\^)}>Be. 

Combining (3.70)-(3.71) with the observation that Xm-i(^m) = V Sg, we 
find that Xm-i{bm) < and hence that hm-i > D 

Proposition 3.11. Suppose that a{-,-) is transient. Then: 

(a) 62 = 2/G(2)(0,0). 

Suppose that a(-,-) is symmetric and transient. Then: 

(b) 62>fe3>64>--->0. 

(c) b2<{m- l)bm < 2/G(™) (0, 0) < 2. 

(d) hm^^oo {m-l)bm=c = sup^>2 (m - 1)6^ • 

Proof, (a) The formula for 62 is obvious because 5^^^ = {0} and Jj^^^ (0) = 
1, giving A2 = G(2)(0,0). 

(b) The fact that m 1-^ is nonincreasing is trivial from (2.1) and 
Lemma 2.4. The fact that bm> for all m follows from (c). 

(c) We prove that 

(3.72) tj("^)(0)G('")(0,0) < A„ < tt(™)(0)G(2)(0,0), 

which implies the claim. The lower bound is obtained by picking n = 60 in 

(2.72) and using (2.75). The upper bound follows from (2.74)-(2.75). 

(d) The proof is an adaptation of the argument in Carmona and Molchanov [7] , 
Chapter III, Section 1.3. The key is the following inequality: 

Lemma 3.12. // 

r 

(3.73) m = ^ Uimi with Ui, mj G N and rrij > 2 for i = 1, . . . ,r, 

1=1 

then 

(3.74) ^J-L-^\<j2na^J^-) V6>0. 
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Proof. Let V denote the collection of all partitions of {!,... ,m} into 
rii groups of mi integers for i = 1, . . . , r. For P G P, write 

(3.75) P= {Pij)i=l,...,r,j=l,...,ni 

to label these groups, so that 



(3.76) 



{l,...,m} = U U^ii 

i=ij=i 



with 

(3.77) \Pij\ = mi, i = l,...,r, j = l,...,ni. 

The cardinality of V is 



(3.78) 

Moreover, 

(3.79) 

with 
(3.80) 
Now, let 
(3.81) 

Then, by (3.73), 
(3.82) 



\V\ = N- 



EE E 1 

lP&Vj=l k,l&P,j 
l<k<l<m 



mi 



Ni, i = l,...,r, 



niiimi — 1) ,^ 
m[m — 1) 



m — 1 1 



mi -1 N 



i=l 



Using (3.79) and (3.82), we may write 

l<k<l<m^ 



(3.83) 



EE'^^E E L ifew=5i(«)} 

P^Vi=l 3=1 k,l&Pij 
l<k<l<m 



ds. 
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nn-P 

\PeVi=i 



. j=i k,i&p,, •'^ 



ds 



l<fc</<m 



(3.^ 



< 



Pev \i=i 



j=i k,ieP„ "'^ J/ 



l<fc</<m 



Pevi=i \ L j=i k,i&Pii J/ 



1/iV 



l<fc</<m 



= n^^'" exp Af65i^ E / lteW=6(«)}'^« ' 
i=i \ L i=i fc,ieP*. ^ \ 

l<fc<Km, 

where P* is any representative partition. Here, the third line uses Holder's 
inequality, the fourth line uses that the factors labelled by i are independent, 
and the fifth line uses that the expectations for given P do not depend on 
the choice of P . Taking logarithms, dividing by t and letting t ^ oo on both 
sides of (3.84), we obtain 



(3.85) 



4 = 1 



[recall (3.40) and (3.77)]. Inserting (3.81) and replacing h by hj{m — 1), we 
arrive at (3.74). □ 

We now show why Lemma 3.12 implies the claim in Proposition 3.11(d). 
An immediate consequence of Lemma 3.12 and the definition of 6m, 

(3.86) 6m = sup{6>0:xm(6) = 0} 

[recall (3.41) and use that hm = bm], is that, subject to (3.73), 



(3.87) 



(m - l)bm > min {mi - 1)6^,- 

1=1,. ..,r 



Let c = limsupm^oo("^ ~ l)6m. Pick e > arbitrary, and pick mi = mi{e) 
such that 



(3.^ 



(mi - l)6mi >c-e. 
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For m > mi, let rii = \m/mi \ . Since m < nimi <m + mi, we may estimate 

Tn — 1 

(m - l)bm > {m - l)6„imi = -{nimi - l)6„imi 

nimi — 1 

(3.89) 

m — 1 , , , m — 1 , , 

> 7 [mi - l)bm, > —— r(c-e). 

nimi — 1 m + mi — 1 

Here, the first inequality uses that m i— > 6^ is nonincr easing, while the second 
inequality uses (3.87) for r = 1. Let m — > oo to obtain that liminfm^oo(m — 
l)^m > c—e. Let e I to get \\m.m^oo{m— l)hm = c. Note that c = sup„>2(m- 
1)6„. □ 



3.5. Proof of Theorem 1.6. (a) Combine (2.1), Proposition 2.5(d) and 
Lemma 3.5. 

(b) Combine (2.18), (3.40), (3.41), Propositions 2.5(c) and 3.8, and the 
fact that Xm{b) is continuous. 

(c) -(d) Combine Propositions 3.10-3.11. 
(e) Use Lemma 3.7 and (3.58). 

4. Survival versus extinction. In Section 4.1 we prove Theorem 1.4. In 
Section 4.2 we argue in favor of Conjecture 1.8. 

4.1. Proof of Theorem 1.4. In Section 4.1.1 we introduce the Palm dis- 
tribution of our process X = {X{t))t>o. In Section 4.1.2 we use the Palm 
distribution to identify 6*. In Section 4.1.3 we prove Theorem 1.4. 

4.1.1. Palm distribution and its stochastic representation. By size-biasing 
our process 

(4.1) X = {X{t))t>o 
at the origin at time T, we obtain a process 

(4.2) X^ = iX^{t))o<t<T 
defined by 

(4.3) dP{{X^{t))o<t<T G •) = ^MT) dP{{X{t))o<t<T G •)• 

In what follows we construct a stochastic representation of in terms of 
a random walk and an SSDE, leading to a stochastic representation of the 
Palm distribution of X, that is, the law of X'^{T). 
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Step 1. Fix r > 0. Let C = (C(^))t>o be the random walk on / with tran- 
sition kernel a(-, •) and jump rate 1, starting at the origin and independent 
of W. Given C, let 

(4.4) xC'^(t) = {iP(t)}.ei 

be the solution of the SSDE 



dXf'' it) = Y: [Xf {t) - X^/ {t)] dt + ^hX^/ {ty dW,{t) 

(4.5) 

+ 6lf '^(t) l{c(T-t)=»} dt, iel,t> 0. 

Here, the difference with the SSDE in (1.2) is the presence of the last term, 
which produces a ^-dependent random potential. Like (1.2), if X^' (0) G £i 
[the space of configurations defined in (1.4)], then (4.5) has a unique strong 
solution on £i with continuous paths [cf. with Theorem 1.1(a)]. Let 

P{-) = law of C on /[°'°°), 

(4-6) 

qC'^(.) = law of X'^'^(T) on [0,oo)^ for given C- 
Then, as we prove in Proposition 4.2 below, 

(4.7) Q^(-) = J Q^'^(-)P((iC) = law of X^(r) on [0,cx))^ 

is the Palm distribution of X. 

Similarly as in Lemma 2.1, we have a representation formula: 

Lemma 4.1. Given a realization of the random walk the process 
(Xf'^(r))T>o starting from X?'^(0) = can be represented as the following 
functional of the Brownian motions: 



rT 

^ J2hiiT-s)=j}dWj{s) 
•^0 je/ 



xf'^(r) = ee-(i/2)^^£;|[^exp 

+ '^mT-s)=aT-s)}ds 

where ^ = {^{t))t>o is the random walk on I with transition kernel a{-, ■) and 
jump rate 1, and the expectation is over ^ conditioned on ^(0) = i and W 
are independent). 



Proof. The proof uses Ito-calculus and is similar to that of Lemma 2.1. 

□ 
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Step 2. The link between X{T) and X'^{T) is given by the following iden- 
tity, which gives an equivalent characterization of the law defined in (4.3): 

Proposition 4.2. For all measurable finite f : [0, ooY M that depend 
on finitely many coordinates, 



(4.9) 



^f'^(/(xC'^(r))) = ii?^(Xo(r)/(x(r))). 



Proof. Write, using Lemma 2.1, 
G 



E^{Xo{T)f{X{T))) 



(4.10) 



: e-(V2)^T^H. L< U r J2 l|c(T-.)=,} dW,{s) 



]f{X{T))] 



exp 



•^0 .e/ 



{C(T-s)=i}d^j(s) 



2 Jo 



T 



l/(^(r)) • 



Define 



(4.11) VF/(t) = P7i(t)-^/6^*l|^(r_,)=i}ds, ieI,t€[0,T]. 

By Girsanov's formula, these are independent standard Brownian motions 
under the transformed path law P^'' defined by 



^pwr(w[o,T]) = exp 



(4.12) 



rT 
■^0 ,v- r 



2 Jo 



J2 l{C(T-s)=j} 



Now, let X^^'^ be the unique solution of the SSDE 
dXP(t) = Y: a(^,i)[x/'^(^) - x/'^(t)] dt 



jG/ 



(4.13) + V^^/'''(i)' dwHt) + &x/'^(t) l{c(T-t)=^} ^^t, 

iel,te [0,T], 

with initial condition X^''^{0) = @. Comparing this with (1.2), we see that 

(4.14) C{X) under = £(1^'^) under P^\ 
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Combining this with (4.10) and (4.12), we obtain 

(4.15) ^E'^{XoiT)f{XiT)))=E<iE'^'ifiX^^'^iT)))). 

The right-hand side of (4.15) is the same as the left-hand side of (4.9). □ 

Step 3. The reason why the Palm distribution is convenient is the following 
important fact, which relies on Proposition 4.2: 

Lemma 4.3. (a) w-limT^oo'C(X^(r)) = Q°^ for some measure Q°° on 
[0,cx)]^. 

(b) Either Q°° = 6oo or Q°° is a probability measure on [0,oo)^. 

(c) w-limT_»oo >C(X(T)) = ^0 o-'^^d only if = 6oo- 

Proof, (a) Following Cox, Klenke and Perkins [11], via the duality re- 
lation in Lemma 2.2 we know that w-limT^oo >C(X(T)) = i'q when X{0) = Q 
[recall (1.22)]. Since E{X(T)) = Q for all T > 0, the Palm measure is defined 
for all T > and converges to a limiting measure Q°° on the compactified 
state space [0,oo]^. Indeed, because X{T) converges in law, for each n G N 
the restriction of X^{T) to An = {x G [0, oo]^ -xq <n} converges in law too. 
The defect of on UneN^ 

n is the weight of oo. 

(b) First we note that the mean of is either G (survival) or (ex- 
tinction). To see why this is true, note the linearity of the SSDE: the so- 
lution for the initial condition QX{0) is {QX(T))t>o- We know from [11], 
Theorem 2.3, that w-limT^oo >C(X(r)) = fe for every shift-invariant and 
shift-ergodic C{X(0)) with mean G. If starting from X(0) = 1 we get a limit 
law with mean a, then by linearity E'^^{Xq) = aQ. However, starting from 
C{X{0)) = vq we then get E'^'^^{Xq) = a^G (since is invariant). Thus, we 
must have = a, that is, a = or 1. 

If vq has mean G, then Xq is uniformly integrable with respect to C{X{T)), 
T > 0, and hence the Palm measures are tight. If, on the other hand, uq, = 5o, 
then Q°° = 5gg_, because all the mass is eventually moving out of An for ev- 
ery n S N. The representation of the Palm measure in (4.8) shows that the 
same holds for all other components as well. 

(c) For a proof we refer to Kallenberg [36]. The intuition is the following. 
First note that, by translation invariance and irreducibility, the statement is 
equivalent to w-limr^oo >C(Xo(T)) = 6o if and only if w-limr^oo >C(X(f (T)) = 
5oo- Next, if X locally dies out, then (because the mean is preserved) it 
clusters, that is, it develops high peaks at rare sites. After size-biasing, as 
in Proposition 4.2, these peaks cause explosion, which is why X diverges. 
Conversely, if X locally survives, then it does not cluster. The size-biasing 
therefore does not cause explosion, which is why X does not diverge. □ 
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From Lemma 4.3 we obtain the following criterion for survival versus 
extinction of X. Local survival means that w-limQ"_»oo>C(X(T)) = vq^ Jq- 

Lemma 4.4. Let P^(-) denote the law of X{T) on [0, cx))^ [recall (4.6) 
and (4.7)]. Then 

{P'^)t>o locally survives <^=^ (Q^)t>o is tight 

kQ )t>o "Is tight C-a.s. 

Proof. The first equivalence in (4.16) uses Lemma 4.3. The second 
equivalence in (4.16) is trivial. □ 

In Section 4.1.2 we use Lemma 4.4 to identify 6*. 

4.1.2. Identification of b^: . Abbreviate 

(4.17) M«(T) = ee"(^/2)feT 

exp 

This is a martingale in W for every fixed ^, with M^(0) = O. In terms of 
this quantity, the representation formulas for X [in (2.2)] and X [in (4.8)] 
read 

x,{t) = eI{mHt)), 

(4.18) 

xP(r) = 4(M5(r)e^^"(^)), 

where T^'''(T) is the collision local time of ^ and C up to time T. 

Proposition 4.5. X survives for < 5 < 6* and locally dies out for 
b> b^, where 

(4.19) 6* = sup{6 > : (g^'^) T>o is tight C-a.s.}. 

Proof. This follows from (4.6) and (4.16). Note that (4.8) shows that 
tightness of (Q^'"^)r>o is measurable with respect to the tail sigma-algebra 
of C, which is trivial. □ 

Equation (4.19) is to be compared with the formulas for bm and 6=k=k given 
in (1.49) and (1.51). However, (4.19) is more difficult to analyze. For one, 
since b appears with both signs in (4.17), it is not a priori obvious that 
b^: defines a unique transition: for this we need to show that if (Q^'"^)t>o 
is tight C-a-S- for some 6 > 0, then the same is true for all smaller values 
of b. Fortunately, the latter property can be shown to hold with the help 
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of a coupling technique put forward in Cox, Fleischmann and Greven [9]. 
There it is shown that, for processes with fixed mean (0 in our case), "more 
noise causes the process to be more spread out and hence to be more prone 
to extinction." More precisely, it is proved that two systems of the type in 
(1.8), with diffusion functions 51, 52 satisfying gi > §2, have the property 

(4.20) ^^.9i(e-^^o(i)) >^x,<;2(g-AXoW) VA>0, 

where E^'^ is expectation over X driven by g. Therefore, if the right-hand 
side converges to 1 as t — > 00, then so does the left-hand side and, conversely, 
if the left-hand side remains bounded away from 1 as t — > 00, then so does 
the right-hand side. Applying this to gi{x) = hix^ and g2{x) = b2x'^, we get 
the required monotonicity. 

Our next observation is the following: 

Lemma 4.6. 6=,, > 6*=,.. 

Proof. Take the expectation over W in the second line of (4.18), use 
Fubini and the fact that (4.17) is a martingale in W, to obtain 

(4.21) £;^(xP(r)) = 4(e^^"(^)). 

Let 6 < 6*=,,. Then, according to (1.51), the right-hand side is C-a.s. bounded 
in T. Consequently, (X^- {T))t>o is a tight family of random variables (^-a.s. 
Thus, by (4.6), (Q^'"^)r>o is tight C-a.s. Therefore, by Lemma 4.4, b< h^. 
□ 

In Section 4.2 we will give an argument suggesting that 6** > 62- To- 
gether with Lemma 4.6 this would imply that 6* > 625 which is the claim in 
Conjecture 1.8. 

4.1.3. Proof of Theorem 1 A. We look at the corresponding parts of The- 
orem 1.3. 

(a)-(b) Combine Proposition 2.5(a) with Proposition 4.5. To get that vq 
has mean 0, see the proof of Lemma 4.3. 

(c)-(d) These follow from Lemma 2.2 along the lines of Cox, Klenke and 
Perkins [11]. The key property is the self-duality expressed in Lemma 2.2. 
Indeed, self-duality translates convergence of the process into convergence 
of the process starting from a configuration with finite mass. For the latter, 
martingale convergence arguments can be applied to get (d). After that, (c) 
follows from (b) and (d). 

(e) This follows from (d) and the representation formula in Lemma 2.1, 
which show that uq tends to 5oo as — > 00 because of the presence of the 
prefactor in (2.2). 
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4.2. Evidence in favor of Conjecture 1.8. In this section we show how 
the conjecture follows from a certain quenched large deviation principle for 
a renewal process in a random environment. This large deviation principle 
is the topic of forthcoming work [5]. 



Step 1. According to (1.50) and (1.51), we have 



)|6> 0:E^'(^exp b J 



(4.22) = sup|0 > U : E'- [ exp 



oo 

bl 1 



dt 



< oo ^-a.s. >. 



Expand the expectation in powers of b to write 



exp 



(4.23) 



N=0 



N 



poo poo 

dsi I ds2 ■ ■ ■ dsN 

J Si •'•5JV-1 



xP«'(e(si)=C'(si),...,e(57v)=e'(^7v)) 



n=0 



with 



(4.24) 



^•00 poo poo 

ds^ I dso ■ ■ dsr, 



Si 



Jo 

X Ps,{0,C{si))Ps,^sMsi),C{s2)) X 
X ^s„-s„-i(e(Sn-l),C(s„)), 



where Pt{i,j) = {S,'{t) = j \ ^'(0) = i), i,j £ I, is the transition kernel of 
the random walk (recall Section 2.5). An easy computation shows that 



(4.25) E^iKniO) = 1 

It therefore follows that 



VnG N. 



(4.26) 



lim -logi^n(C) = C< C-a.s., 

n~*oo n 



where the fact that the limit is .^-a.s. constant is an immediate consequence 
of the fact that the exponential growth rate of K^i^,) is measurable with 
respect to the tail sigma-algebra of ^, which is trivial. By (4.22)-(4.23), to 
get > &2 it suffices to show that C < 0. 



BRANCHING PROCESS IN A RANDOM ENVIRONMENT 55 

Step 2. Let {Tk)k&'} be i.i.d. positive random variables with probability 
law 

Since 62 = 2/G(0,0) by (1.20), this gives 
where 

(4.29) Afc = e[^r„Ve[^r„Y A; G N, 

\m,=l / \m=l / 

with tq = 0. The expectation in (4.28) seems well suited for a large deviation 
analysis, but the problem is that the are correlated because ^ is fixed. 

Step 3. Abbreviate 

k 

(4.30) n = Y.ru n = C[T,_„T,), ken. 

1=1 

Then {Yk)keN is a random process taking values in Q = {J^yQD{[0,t),I), 
with D{[0,t),I) the Skorohod space of paths of length t. For z — ^[o,t) 
we abbreviate 

(4.31) Az = z{t) - z{0), \z\=t-0 = t. 
Introduce the empirical process 



1 

(4.32) Rn = -Yl kv. 



" fc = l 



n G 



where the upper index per means that the n-sequence is periodically re- 
peated. For each G N, we have G O. Consequently, Rn is a random ele- 
ment of V{^)^ the set of probability measures on O^, with the randomness 
coming from (rfc)fcgN (C being fixed). In terms of Rn we may write 



(4.33) i^„(0=^"^'-'"" exp 



n / {■KiRn){dy)\or 



In' " ^^2^(0,0) 

where iriRn is the projection of i?„ onto the first coordinate. 
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Step 4. In [5] it is shown that (-Rn)neN satisfies a quenched (with respect 
to ^) large deviation principle on 7^(Q^) with some rate function Q i— > I{Q). 
This rate function turns out to be a sum of two terms, the first being the 
specific relative entropy of Q with respect to Q^, where 

(4.34) QO = £((n,)feGN), 

the second being a specific relative entropy associated with "concatenation 
and randomization" of (yfc)fcgN. By applying Varadhan's lemma to (4.33), 
we therefore obtain 

(4.35) lira -log Kn(0=C, ^a.s. 
with 



(4.36) C = sup 



/KQ)(^^)iog^^-/(g) 

Jn F2\y\[0,u) 



We note that Varadhan's lemma applies because {t,i) i— > logPf(0,i)/-P2i(0,0) 
is bounded from above and tending to — oo as ||i|| — > oo (see Deuschel and 
Stroock [21], Theorem 2.1.10). 

Step 5. Note that is a product measure: 

(4.37) = (gO)^ with q^idz) = q\dT, ^^[0,^)) = i^(dr)/.(dC[o,r))- 
Since the integral in (4.36) depends on vriQ only, by putting 

(4.38) inf I{Q)=i{q), 

niQ=q 

we can reduce the variational expression in (4.36) to 



(4.39) C= sup 



l{dy)log ' ' -i{q) 



This formula is the key to proving that C < 0. Indeed, i{q) turns out to be a 
sum of two terms, one being h{q\q^) (the relative entropy of q with respect 
to q^), the other being a relative entropy term associated with concatenation 
and randomization. Without the second term, the supremum in (4.39) would 
be attained at 

(4.40) q{dz) = q{dT, ^^[o,,)) = 62^,(6, ^(r)) dr fi{d^[o,r)) 

and would be equal to 0. Indeed, this corresponds to the annealed upper 
bound in (4.25) and (4.26). With the second term, however, it is < 0. See 
[5] for further details. 
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